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Abstract: Expectation values of Wilson loops define the nonperturbative properties of the hot 
medium produced in heavy ion collisions that arise in the analysis of both radiative parton energy 
loss and quarkonium suppression. We use the AdS/CFT correspondence to calculate the expectation 
values of such Wilson loops in the strongly coupled plasma of Af = 4 super Yang-Mills (SYM) theory, 
allowing for the possibility that the plasma may be moving with some collective flow velocity as is 
the case in heavy ion collisions. We obtain the M = 4 SYM values of the jet quenching parameter 
q, which describes the energy loss of a hard parton in QCD, and of the velocity-dependence of 
the quark-antiquark screening length for a moving dipole as a function of the angle between its 
velocity and its orientation. We show that if the quark-gluon plasma is flowing with velocity Vf 
* at an angle 9 with respect to the trajectory of a hard parton, the jet quenching parameter q is 
modified by a factor 7^ (1 — Vf cos#), and show that this result applies in QCD as in M = 4 SYM. 
We discuss the relevance of the lessons we are learning from all these calculations to heavy ion 
collisions at RHIC and at the LHC. Furthermore, we discuss the relation between our results and 
those obtained in other theories with gravity duals, showing in particular that the ratio between q in 
any two conformal theories with gravity duals is the square root of the ratio of their central charges. 
This leads us to conjecture that in nonconformal theories q defines a quantity that always decreases 
along renormalization group trajectories and allows us to use our calculation of q in M = 4 SYM 
to make a conjecture for its value in QCD. 
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1. Introduction 



Understanding the implications of data from the Relativistic Heavy Ion Collider (RHIC) poses 
qualitatively new challenges [l|]. The characteristic features of the matter produced at RHIC, namely 
its large and anisotropic collective flow and its strong interaction with (in fact not so) penetrating 
hard probes, indicate that the hot matter produced in RHIC collisions must be described by QCD 
in a regime of strong, and hence nonperturbative, interactions. In this regime, lattice QCD has to 
date been the prime calculational tool based solely on first principles. On the other hand, analyzing 
the very same RHIC data on collective flow, jet quenching and other hard probes requires real-time 
dynamics: the hot fluid produced in heavy ion collisions is exploding rather than static, and jet 
quenching by definition concerns probes of this fluid which, at least initially, are moving through it 
at close to the speed of light. Information on real-time dynamics in a strongly interacting quark- 
gluon plasma from lattice QCD is at present both scarce and indirect. Complementary methods 
for real-time strong coupling calculations at finite temperature are therefore desirable. 

For a class of non-abelian thermal gauge field theories, the AdS/CFT conjecture provides such 
an alternative f2|. It gives analytic access to the strong coupling regime of finite temperature gauge 
field theories in the limit of large number of colors (N c ) by mapping nonperturbative problems at 
strong coupling onto calculable problems in the supergravity limit of a dual string theory, with the 
background metric describing a curved five- dimensional anti-deSitter spacetime containing a black 
hole whose horizon is displaced away from "our" 3+1 dimensional world in the fifth dimension. 
Information about real-time dynamics within a thermal background can be obtained in this set- 
up. The best-known example is the calculation of the shear viscosity in several supersymmetric 
gauge theories || f|, ^ [7], ^ |], |TU|. It was found that the dimensionless ratio of the shear 
viscosity to the entropy density takes on the "universal" @, [|, ||, |10| value 1/4tt in the large 
number of colors (N c ) and large 't Hooft coupling (A = g\ M N c ) limit of any gauge theory that 
admits a holographically dual supergravity description. Although the AdS/CFT correspondence 
is not directly applicable to QCD, the universality of the result for the shear viscosity and its 
numerical coincidence with estimates of the same quantity in QCD made by comparing RHIC data 
with hydrodynamical model analyses |Tl| have motivated further effort in applying the AdS/CFT 
conjecture to calculate other quantities which are of interest for the RHIC heavy ion program. This 
has lead to the calculation of certain diffusion constants 0] and thermal spectral functions as 
well as to first work (|I4J towards a dual description of dynamics in heavy ion collisions themselves. 
More recently, there has been much interest in the AdS/CFT calculation of the jet quenching 
parameter which controls the description of medium-induced energy loss for relativistic partons in 
QCD [16|, [17|, [18], [19], [20|, ^ ^2) and the drag coefficient which describes the energy loss for 

There have also 



heavy quarks in M = 4 supersymmetric Yang-Mills theory [p3 , 24[ ^5 , |26| , |27 



been studies of the stability of heavy quark bound states in a thermal environment [28, [29], [30]] with 
collective motion [BO, Bl|, [ 



33, 34, 35, pa m, 38, 39 



The expectation values of Wilson loops contain gauge invariant information about the nonper- 
turbative physics of non-abelian gauge field theories. When evaluated at temperatures above the 
crossover from hadronic matter to the strongly interacting quark-gluon plasma, they can be related 
to a number of different quantities which are in turn accessible in heavy ion collision experiments. 
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In Section 2 of this paper, which should be seen as an extended introduction, we review these con- 
nections. We review how the expectation value of a particular time-like Wilson loop, proportional 
to exp(— iS) for some real S, serves to define the potential between a static quark and antiquark 
in a (perhaps moving) quark-gluon plasma. However, in order to obtain a sensible description of 
the photo-absorption cross-section in deep inelastic scattering, the Cronin effect in proton-nucleus 
collisions, and radiative parton energy loss and hence jet quenching in nucleus-nucleus collisions, 
the expectation value of this Wilson loop must be proportional to exp(— S) for some real and 
positive S once the Wilson loop is taken to lie along the lightcone. In Section 3, we present the 
calculation of the relevant Wilson loops in hot M = 4 supersymmetric Yang- Mills theory, using the 
AdS/CFT correspondence, and show how its expectation value goes from exp(— iS) to exp(— S) (as 
it must if this theory is viable as a model for the quark-gluon plasma in QCD) as the order of two 
non-commuting limits is exchanged. The jet quenching parameter q, which describes the energy 
loss of a hard parton in QCD, and the velocity-dependent quark-antiquark potential for a dipole 
moving through the quark-gluon plasma arise in different limits of the same Nambu-Goto action 
which depends on the dipole rapidity rj and on A, the location in r, the fifth dimension of the AdS 
space, of the boundary of the AdS space where the dipole is located. If we take r] — > oo first, and 
only then take A — ► oo, the Nambu-Goto action describes a space-like world sheet bounded by a 
light-like Wilson loop at r = A, and defines the jet quenching parameter. If instead we take A — > oo 
first, the action describes a time-like world sheet bounded by a time-like Wilson loop, and defines 
the gg-potential for a dipole moving with rapidity 77. We review the calculation of both quantities. 
In Section 4 we calculate the jet quenching parameter in a moving quark-gluon plasma, and show 
that our result in this section is valid in QCD as in M = 4 SYM. In Section 5 we return to the 
velocity-dependent screening length, calculating it for all values of the angle between the velocity 
and orientation of the quark-antiquark dipole. 

Section 6 consists of an extended discussion. We summarize our results on the velocity- 
dependent screening length in Section 6.1. In Section 6.2, we comment on the differences between the 
calculation of the jet quenching parameter and the drag force on a (heavy) quark |23|, [24], [2l| 27f . 



In Section 6.3, we then compare our calculation of the jet quenching parameter to the value of this 
quantity extracted in comparison with RHIC data. The success of this comparison motivates us 
to, in Section 6.4, enumerate the differences between QCD and M = 4 supersymmetric Yang-Mills 
(SYM) theory, which have qualitatively distinct vacuum properties, and the rapidly growing list of 
similarities between the properties of the quark-gluon plasmas in these two theories. A comparison 
between our result for the velocity scaling of the quark-antiquark screening length and future data 
from RHIC and the LHC on the suppression of high transverse momentum J/ty and T mesons 
could add one more entry to this list. The single difference between M = 4 SYM and QCD which 
appears to us most likely to affect the value of the jet quenching parameter is the difference in the 
number of degrees of freedom in the two theories. We therefore close in Section 6.5 by reviewing the 
AdS / CFT calculations to date of q in theories other than M = 4 SYM, and show that for any two 
conformal field theories in which this calculation can be done, the ratio of q in one theory to that in 
the other will be given by the square root of the ratio of the central charges, and hence the number 
of degrees of freedom. This suggests that q in QCD is smaller than that in M = 4 SYM by a factor 
of order a/120/47.5 ~ 1.6. This conjecture can be tested by further calculations in nonconformal 
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Figure 1: Schematic illustration of the shape of Wilson loops C, corresponding to a qq dipole of size 
L, oriented along the xi-direction, which is (i) at rest with respect to the medium (C sta tic ) , (h) moving 
with some finite velocity v = tanhry along the longitudinal ^-direction (Ck t °°?* ed ), or ( m ) moving with the 
velocity of light along the ^-direction (Ciight-iike)- 

theories. 

A reader interested in our results and our perspective on our results should focus on Sections 
4 and 6. A reader interested in how we obtain our results should focus on Section 3. 



2. Wilson loops in heavy ion collisions 

In this section, we consider Wilson lines 



W r {C) = TrPexp 



h I dx fl A fl (x) 



(2.1) 



where J c denotes a line integral along the closed path C. W r (C) is the trace of an SU(N)-mairix 
in the fundamental or adjoint representation, r = F, A, respectively. The vector potential A^(x) = 
A%(x) T a can be expressed in terms of the generators T a of the corresponding representation, and 
V denotes path ordering. We discuss several cases in which nonperturbative properties of interest 
in heavy ion physics and high energy QCD can be expressed in terms of expectation values of (|2.1[ ). 



2.1 The quark-antiquark static potential 

We shall use the Wilson loop 

(W F (C Btatix >)) = exp [-iT (E(L) - E rcn )} 



(2.2) 



to furnish a working definition of the qq static potential E(L) for an infinitely heavy quark-antiquark 
pair at rest with respect to the medium and separated by a distance L. Here, the closed contour 
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Static has a short segment of length L in the transverse direction, and a very long extension T in the 
temporal direction, see Fig. [I]. The potential E(L) is defined in the limit T — > oo. The properties 
of the medium, including for example its temperature T, enter into (|2.2|) via the expectation value 
(...). Here, E ren is an L-independent renormalization, which is typically infinite. Eq. (|2.2| ) is 
written for a Minkowski metric, as appropriate for our consideration, below, of a quark- ant iquark 
pair moving through the medium. At zero temperature, the analytic continuation iT — > T of (|2.2| ) 
yields the standard relation between the static potential and an Euclidean Wilson loop jl(J. In 
finite temperature lattice QCD J|I], 0, f43fl , one typically defines a quark- ant iquark static potential 



from the correlation function of a pair of Polyakov loops wrapped around the periodic Euclidean 
time direction. (For a discussion of this procedure and alternatives to it, see also Ref. f44|.) In 
these Euclidean finite temperature lattice calculations, the corresponding quark-antiquark potential 
is renormalized such that it matches the zero temperature result at small distances [fi3| |. We shall 



use an analogous prescription. We note that while (|2.2|) is difficult to analyze in QCD, its evaluation 
is straightforward for a class of strongly interacting gauge theories in the large number of colors 



limit at both zero [45| and nonzero temperature [ 28| , as we shall see in Section 3. 



The dissociation of charmonium and bottomonium bound states has been proposed as a sig- 
nal for the formation of a hot and deconfined quark-gluon plasma [Q. Recent analyses of this 
phenomenon are based on the study of the quark-antiquark static potential extracted from lattice 
QCD |I7||. In these calculations of E(L), the gg-dipole is taken to be at rest in the thermal medium, 
and its temperature dependence is studied in detail. In heavy ion collisions, however, quarkonium 
bound states are produced moving with some velocity v = tanhr? with respect to the medium. If 
the relative velocity of the quarkonium exceeds a typical thermal velocity, one may expect that 
quarkonium suppression is enhanced compared to thermal dissociation in a heat bath at rest [pT 



For a calculation of the velocity-dependent dissociation of such a moving gg-pair in a medium at 
rest in the a^-direction, one has to evaluate Q2.2|) for the Wilson loop C^°^ ed , depicted in Fig. |l[ 
The orientation of the loop in the (x^, t)-plane changes as a function of rj. This case is discussed in 
section |3.1| . In section |5|, we discuss the generalization to dipoles oriented in an arbitrary direction 
in the (xi, X3)-plane. 



2.2 Eikonal propagation 

We now recall cases of physical interest where, unlike in (|2.2| ), the expectation value of a Wilson 
loop in Minkowski space is the exponent of a real quantity. Such cases are important in the high 
energy limit of various scattering problems. Straight light-like Wilson lines of the form H^Xj) = 
V exp{i J dz~T a A'^('Xi, z~)} typically arise in such calculations when — due to Lorentz contraction 
- the transverse position of a colored projectile does not change while propagating through the 
target. The interaction of the projectile wave function with the target can then be described in 
the eikonal approximation as a color rotation aj — > fa of each projectile component i, resulting in 
an eikonal phase W ai/ 3 i ('Xi) . A general discussion of this eikonal propagation approximation can 
be found in Refs. Here, we describe two specific cases, in which expectation values of a 

fundamental and of an adjoint Wilson loop arise, respectively. 
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2.2.1 Virtual photoabsorption cross section 

In deep inelastic scattering (DIS), a virtual photon 7* interacts with a hadronic target. At small 
Bjorken x, DIS can be formulated by starting from the decomposition of the virtual photon into 
hadronic Fock states and propagating these Fock states in the eikonal approximation through the 
target J5D], |5l], |52], |53], |54| . However, in a DIS scattering experiment the virtual photon does not have 
time to branch into Fock states containing many soft particles (equivalently, it does not have time 
to develop a colored field) prior to interaction, as it would if it could propagate forever. Instead, the 
dominant component of its wave function which interacts with the target is its qq Fock component: 



1 * \ 
17 ) 



d 2 (x - y) dzip(x - y, z)— =S a& \a(yi) , a(y), z) . 



N 



(2.3) 



Here, |a(x) , a(y), z) denotes a gg-state, where a quark of color a carries an energy fraction z 
and propagates at transverse position x. The corresponding antiquark propagates at transverse 
position y and carries the remaining energy. The Kronecker 5 aSl ensures that this state is in a color 
singlet. A" is the number of colors; the probability that the photon splits into a quark antiquark 
pair with any one particular color is proportional to 1/N. The wave function ip is written in the 
mixed representation, using configuration space in the transverse direction and momentum space 
in the longitudinal direction. It can be calculated perturbatively from the 7* — > qq splitting |53 . 
Given an incoming state \^- m ) = l a ( x ) ; «(y)), in the eikonal approximation the outgoing state 
reads l^out) = W^( x ) ^J=/(y)|7( x ) , 7(y))> an d the total cross section is obtained by squaring 
l^tot) = l^out) — l^m)- From the virtual photon state Q2.3|), one finds in this way the total virtual 



photoabsorption cross section [58 



a 



DIS 



p39 
-'tot 



d 2 x d 2 y dzip(pc — y , z) tp* (x - y , z) P™ t (x, y ) , 
iTr [W f {x) W F ^(y)} - iTr [W F (y) W F \x 



N 



(2.4) 



(2.5) 



This DIS total cross section is written in terms of the expectation value of a fundamental Wilson 
loop: 



1 

N 



(Tr [W F \y) W F (x)}) — <W F (C ligh t-i ik e)> = exp 



1 



Q S L 



+ 



1 

AT2 



(2.6) 



By the — > we mean that in order to obtain a gauge-invariant formulation, we have connected the 
two long light-like Wilson lines separated by the small transverse separation L = |x — y| by two 
short transverse segments of length L, located a long distance L~ ^> L apart. This yields the closed 
rectangular loop Ci ig ht-iikc illustrated in Fig. |l|. The expectation value (...) denotes an average over 
the states of the hadronic target; technically, this amounts to an average over the target color fields 
in the Wilson line ( [2.1|) . If we could do deep inelastic scattering off a droplet of quark-gluon 
plasma, the (...) would be a thermal expectation value. We have parameterized (H^ F (Cii g ht-iike)) in 
terms of the saturation scale Q 2 S . This is the standard parametrization of virtual photoabsorption 
cross sections in the saturation physics approach to DIS off hadrons and nuclei [pR |~~ 



57]. Although 
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we do not know the form of the 1/N 2 corrections to (|2.6|) , we do know that they must be such that 
{W F } -> in the L -> limit, since in this limit P^ t must vanish. 

We note that for small L, the /^-dependence of the exponent in ( |2.(j| ) follows from general 
considerations. Since the transverse size of the gg-dipole is conjugate to the virtuality Q of the 
photon, L 2 ~ 1/Q 2 , one finds P™ t = \Q 2 S L 2 + 0(L 4 ) ~ Q 2 JQ 2 . This is the expected leading 
Q 2 -dependence at high virtuality. 



General considerations also indicate that the exponent in (2.6) must have a real part. To see 
this, consider the limit of large L and small virtuality, when the virtual photon is large in transverse 
space, and its local interaction probability should go to unity. Since Eq. (|2.5|) is the sum of the 



elastic and inelastic scattering probability, which are both normalized to one, one requires P^ t — > 2 
in this large- L limit. This cannot be achieved with an imaginary exponent in ( |2.6| ). 

The saturation momentum Q s is a characteristic property of any hadronic target. Qualitatively, 
the gluon distribution inside the hadronic target is dense (saturated) as seen by virtual photons up 
to a virtuality Q s , but it is dilute as seen at higher virtuality. As a consequence, a virtual photon 
has a probability of order one for interacting with the target, if — in a configuration space picture 
- its transverse size is |x — y| > 1/Q S , and it has a much smaller probability of interaction for 



y| <C l/Qs- This is the physics behind ( |2.4| ) and 



2.2.2 The Cronin effect in proton-nucleus (p-A) collisions 

In comparing transverse momentum spectra from proton-nucleus and proton-proton collisions, one 
finds that in an intermediate transverse momentum range of pr ~ 1 — 5 GeV, the hadronic yield 



in p-A collisions is enhanced |p~8|| . This so-called Cronin effect is typically understood in terms of 
the transverse momentum broadening of the incoming partons in the proton projectile, prior to 
undergoing the hard interaction in which the high-p^ parton is produced. On the partonic level, 
this phenomenon and its energy dependence have been studied by calculating the gluon radiation 
induced by a single quark in the incoming proton projectile scattering on a target of nuclear size A 
and corresponding saturation scale Q s [^9], [61], [62], |63|, |64| • 

One starts from the incoming wave function ty? n of a bare quark |a(0)), supplemented by the 
coherent state of quasi- real gluons which build up its Weizsacker- Williams field /(x) cx g 2 ^. Here, 
g is the strong coupling constant and x and are the transverse positions of the gluon and parent 
quark fjgfl . Suppressing Lorentz and spin indices, one has = \a(0))+j (htd^ /(x) T^|/5(0); 6(x, £)). 
The ket |/3(0); 6(x, £)) describes the two-parton state, consisting of a quark with color (3 at trans- 
verse position and a gluon of color b at transverse position x. In the eikonal approximation, the 
distribution of the radiated gluon is flat in rapidity £. The outgoing wave function differs from tyf n 
by color rotation with the phases W£p for quarks and W^ c for the gluons: 

*out = W£(0)| 7 > + J rfx/(x)T^Mf 7 (0)^ c (x)| 7 ;c(x)). (2.7) 

(a, (3 and 7 are fundamental indices; b, c and d below are adjoint indices.) To calculate an observable 
related to an inelastic process, such as the number of gluons dN pro< i/dk produced in the scattering, 
one first determines the component of the outgoing wave function, which belongs to the subspace 
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orthogonal to the incoming state \5^f) = [1 — ^^{^nj] \^ ou t)- Next, one counts the number of 
gluons in this state |J5], (KJ 



a.d 



~dk N 



2tt J y x 2 y 2 iV 2 -l 



(Tr [W A \0) W A (0)] ) - <Tr [W a \k) W a {0)] > 



(Tr [W A \y) W A {0)} ) + (Tr [W A \y) W A (x)} ) 



(2.8) 



Here, x and y denote the transverse positions of the gluon in the amplitude and complex conjugate 
amplitude. The fundamental Wilson lines W F (0) at transverse position 0, which appear in (|2.7|) , 
combine into an adjoint Wilson line via the identity W A b (0) = 2Tr \W F {Q>)T a W F] {0)T b ] . We now 
see that the only information about the target which enters in (|2.8|) is that encoded in the transverse 
size dependence of the expectation value of two light-like adjoint Wilson lines, which we can again 
close to form a loop: 



1 



N 2 - 1 



(Tr [W A \y) W A {*)] ) — ► (W A (C light -x^)) = exp 



+ o 



1 

iV 2 



(2.9) 



the parameterization of the expectation 
and fl2.9|) respectively differs in the 



Consistent with the identity Tr W A = Tr W F Tr W F - 1 
values of the adjoint and fundamental Wilson loops in 
large- N limit only by a factor of 2 in the exponent. 

Inserting ( |2.9| ) into ( |2.8| ), Fourier transforming the Weizsacker- Williams factors and doing the 
integrals, one finds formally 



dK 



prod 



dk 



47T 

Q 



— I dq exp 



q- 



Ql\ k 2 (q-k)^ 



(2.10) 



To interpret this expression, we recall the high energy limit for gluon radiation in single quark-quark 
scattering. For a transverse momentum transfer q between the scattering partners, the spectrum in 

2 

the gluon transverse momentum k is proportional to the so-called Bertsch-Gunion factor k 2(q_ k ^ • 
Hence, Eq. ( |2.10p indicates that the saturation scale Q s characterizes the average squared transverse 
momentum q 2 transferred from the hadronic target to the highly energetic partonic projectile. We 
caution the reader that the integrals in (|2.8j ) are divergent and that the steps leading to ( p,10|) remain 
formal since they were performed without proper regularization of these integrals. Furthermore, a 
more refined parametrization of the saturation scale in QCD includes a logarithmic dependence of Q s 
on the transverse separation L. Including this correction allows for a proper regularization |65|, Q . 
The analysis of (|2.8j ) is more complicated, but the lesson drawn from ( [2.10 ) remains unchanged: 
the saturation scale Q 2 S determines the average squared transverse momentum, transferred from the 
medium to the projectile. 
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The dependence of the saturation scale Q 2 on nuclear size A is Q 2 oc A 1 / 3 , i.e., Q 2 is linear 
in the in-medium path- length. According to (|2.1Up , transverse momentum is accumulated in the 
hadronic target due to Brownian motion, q 2 oc A 1 / 3 . In the discussion of high-energy scattering 
problems in heavy ion physics, where the in-medium path length depends on the geometry and 
collective dynamics of the collision region, it has proven advantageous to separate this path-length 
dependence explicitly |)6 

Q" = ^ (2- 11 ) 

in so doing defining a new parameter q. Here, we have expressed the longitudinal distance Az = 
in terms of the light-cone distance L~ . The parameter q characterizes the average transverse 
momentum squared transferred from the target to the projectile per unit longitudinal distance 
travelled, i.e. per unit path length. Note that q is well-defined for arbitrarily large L~ in an infinite 
medium, whereas Q 2 S diverges linearly with L~ and so is appropriate only for a finite system. We 



shall see in Section 2.3 that, when the expectation value in (2^9) is evaluated in a hot quark-gluon 
plasma rather than over the gluonic states of a cold nucleus as above, the quantity q governs the 
energy loss of relativistic partons moving through the quark-gluon plasma. The simpler examples 
we have introduced here in Section 2.2 motivate the need for a nonperturbative evaluation of the 
light-like Wilson loop (^(Cnght-nke)) in a background corresponding to a hadron or a cold nucleus, 
as in so doing one could calculate the saturation scale and describe DIS at small x and the Cronin 
effect. Unfortunately, although hot M = 4 supersymmetric Yang- Mills theory describes a system 
with many similarities to the quark-gluon plasma in QCD as we shall discuss in Section 6, it does 
not seem suited to modelling a cold nucleus. 

2.3 BDMPS radiative parton energy loss and the jet quenching parameter 

In the absence of a medium, a highly energetic parton produced in a hard process decreases its virtu- 
ality by multiple parton splitting prior to hadronization. In a heavy ion collision, this perturbative 
parton shower interferes with additional medium-induced radiation. The resulting interference pat- 



tern resolves longitudinal distances in the target |6§, |69fl . As a consequence, its description 
goes beyond the eikonal approximation, in which the entire target acts totally coherently as a single 
scattering center. As we shall explain now, this refined kinematical description does not involve ad- 
ditional information about the medium beyond that already encoded in the jet quenching parameter 
q that we have already introduced. 



In the Baier-Dokshitzer-Mueller-Peigne-Schiff |67|] calculation of medium-induced gluon radia 



tion, the radiation amplitude for the medium- modified splitting processes q ^ q g 01 g ^ g g is 
calculated for the kinematic region 

£?>w»|k|,|q| = |^q i |>r,A QCD . (2.12) 

i 

The energy E of the initial hard parton is much larger than the energy u of the radiated gluon, which 
is much larger than the transverse momentum k of the radiated gluon or the transverse momentum 
q accumulated due to many scatterings of the projectile inside the target. This ordering is also 
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at the basis of the eikonal approximation. In the BDMPS formalism, however, terms which are 
subleading in 0(1/ E) are kept and this allows for a calculation of interference effects. To keep 
0(l/i?)-corrections to the phase of scattering amplitudes, one replaces eikonal Wilson lines by the 
retarded Green's functions IBHI FFOl |6£l 



G{x 2 ,r 2 ; x x ,ri|p) 



r(x 2 )=r 2 



Dr(x ) exp 



r(x 1 )=n 



dx 



dr(x 
dx~ 



— i 



dx A + (x ,r(x )) 



(2.13) 

Here, p is the total momentum of the propagating parton, and the color field A + = T a is in the 
representation of the parton. The integration goes over all possible paths r(x~) in the light-like 
direction between ri = r(x~[) and r 2 = rix"^)- Green's functions of the form (|2.13|) are solutions to 



the Dirac equation in the spatially extended target color field A + |71| , |70| , p4| . In the limit of ultra- 
relativistic momentum p — > oo, Eq. ( |2.13| ) reduces to a Wilson line ( |2.1| ) along an eikonal light-like 
direction. In the BDMPS formalism, the inclusive energy distribution of gluon radiation from a 
high energy parton produced within a medium can be written in terms of in-medium expectation 
values of pairs of Green's functions of the form ( |2.13| ), one coming from the amplitude and the other 
coming from the conjugate amplitude. After a lengthy but purely technical calculation, it can be 

m 



written in the form 
dl 



LO- 



R 



duj dk 



(2tt) 2 lu 2 

d d 
dy du 



2Re / dyi / dyi due 



pu=r (yi) 




/ Vr exp 




/y=0 


i J yi x 



-iku 



exp 



my 



(2.14) 



Here, the Casimir operator Cr is in the representation of the parent parton. In the configuration 
space representation used in ( |2.14|) , £o is the position at which the initial parton is produced in 
a hard process and the internal integration variables yi and y~\ denote the longitudinal position at 
which this initial parton radiates the gluon in the amplitude and complex conjugate amplitude, 
respectively. (See Refs. pUI, EDI for details.) Since all partons propagate with the velocity of light, 
these longitudinal positions correspond to emission times yi, y"\. 
In deriving ( |2.14|) 



69 



9J, the initial formulation of the q — ► qg radiation amplitude of course 
involves Green's functions ( |2.13j ) in both the fundamental and in the adjoint representation. How- 
ever, via essentially the same color algebraic identities which allowed us to write the gluon spectrum 
(TE) in terms of expectation values of adjoint Wilson loops only, the result given in (|2.14|) has been 



written in terms of expectation values of adjoint light-like Green's functions of the form ( |2.13|) only. 
These in turn have been written in terms of the same jet quenching parameter q defined as in ( |2.9j ) 
and ( |2.11| ), namely via [f49 



(^(Ciight-uke)) = exp 



qL~ L 2 



4V2 



(2.15) 



now with the expectation value of the light-like Wilson loop evaluated in a thermal quark-gluon 



plasma rather than in a cold nucleus. The quantity q(£) which arises in ( 2.14 ) is the value of q at 
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the longitudinal position £, which changes with increasing £ as the plasma expands and dilutes. In 
our analysis of a static medium, = q is constant. 

In QCD, radiative parton energy loss is the dominant energy loss mechanism in the limit in which 
the initial parton has arbitrarily high energy. To see this, we proceed as follows. Note first that in 
this high parton energy limit the assumptions ( |2.12| ) underpinning the BDMPS calculation become 
controlled. And, given the ordering of energy scales in Q2.12|) , the quark-gluon radiation vertex 
should be evaluated with coupling constant a s (k. 2 ). The distribution of the transverse momenta 
of the radiated gluon is peaked around k 2 ~ Q 2 S = qL~ / \/2 Jf2| which means that, in the limit of 
large in- medium path length L~ / y/2, the coupling a s is evaluated at a scale k 2 ^> T 2 at which it is 
weak, justifying the perturbative BDMPS formulation |)7]. Next, we note that in the limit of large 
in-medium path length the result Q2.14 ) yields [67, 73 



dl olJ3 
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R 



du 



2 Re In 



7T 



COS 



(2.16) 



Integrating this expression over lu, one finds that the average medium-induced parton energy loss 
is given by 

1 L- 2 

AE = -a s C R q— , (2.17) 

which is independent of E and quadratic in the path length L~ } This makes the energy lost by 
gluon radiation parametrically larger in the high energy limit than that lost due to collisions alone, 
which grows only linearly with path length, and makes radiative energy loss dominant in the high 
parton energy limit. Radiative parton energy loss has been argued to be the dominant mechanism 
behind jet quenching at RHIC [|68| , |69| , |75| , [76]| , where the high energy partons whose energy loss 
is observed in the data have transverse momenta of at most about 20 GeV p|. At the LHC, the 
BDMPS calculation will be under better control since the high energy partons used to probe the 
quark-gluon plasma will then have transverse momenta greater than 100 GeV [fnj . 

Although the BDMPS calculation itself is under control in the high parton energy limit, a 
weak coupling calculation of the jet quenching parameter q is not, as we now explain. Recall that 
q is the transverse momentum squared transferred from the medium to either the initial parton 
or the radiated gluon, per distance travelled. In a weakly coupled quark-gluon plasma, in which 
scatterings are rare, q is given by the momentum squared transferred in a single collision divided 
by the mean free path between collisions. Even though the total momentum transferred from the 
medium to the initial parton and to the radiated gluon is perturbatively large since it grows linearly 
with the path length, the momentum transferred per individual scattering is only of order g(T) T. 
So, a weak-coupling calculation of q is justified only if T is so large that physics at the scale T is 
perturbative. Up to a logarithm, such a weak-coupling calculation yields p7|, |66|, 



Qweak— coupling 



*C(3) 



-odN T 



71 



(2.18) 



1 For any finite L , corrections to ( 2.14 ) can make the average energy loss AE grow logarithmically with E at 
large enough E [fT4|. 
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if N, the number of colors, is large. However, given the evidence from RHIC data [Tp (the mag- 
nitude of jet quenching itself; azimuthal anisotropy comparable to that predicted by zero-viscosity 
hydrodynamics) that the quark-gluon plasma is strongly interacting at the temperatures accessed 
in RHIC collisions, there is strong motivation to calculate q directly from its definition via the light- 
like Wilson loop ( |2.15| ), without assuming weak coupling. If and when the quark-gluon plasma is 
strongly interacting, the coupling constant involved in the multiple soft gluon exchanges described 
by the weak-coupling calculation of q is in fact nonperturbatively large, invalidating fl2.18| ). 

To summarize, the BDMPS analysis of a parton losing energy as it traverses a strongly inter- 
acting quark-gluon plasma is under control in the high parton energy limit, with gluon radiation 
the dominant energy loss mechanism and the basic calculation correctly treated as perturbative. In 
this limit, application of strong coupling techniques to the entire radiation process described by Eq. 
fl2.14| ) would be inappropriate, because QCD is asymptotically free. The physics of the strongly 
interacting medium itself enters the calculation through the single jet quenching parameter q, the 
amount of transverse momentum squared picked up per distance travelled by both the initial parton 
and the radiated gluon. A perturbative calculation of q is not under control, making it worthwhile 
to investigate any strong coupling techniques available for the evaluation of this one nonperturbative 
quantity. 



3. Wilson loops from AdS/CFT in J\f = 4 super Yang-Mills theory 

In section ^|, we have recalled measurements of interest in heavy ion collisions, whose description 
depends on thermal expectation values of Wilson loops. For questions related to the dissociation 
of quarkonium, the relevant Wilson loop is time-like and (W(C)) is the exponent of an imaginary 
quantity. Questions related to medium-induced energy loss involve light-like Wilson loops and 
(W(C)) is the exponent of a real quantity. 

In this section, we evaluate thermal expectation values of these Wilson loops for thermal M = 
4 super Yang-Mills (SYM) theory with gauge group SU(N) in the large N and large 't Hooft 
coupling limits, making use of the AdS/CFT correspondence [0, In the present context, this 



correspondence maps the evaluation of a Wilson loop in a hot strongly interacting gauge theory 
plasma onto the much simpler problem of finding the extremal area of a classical string world sheet 
in a black hole background p8| |. We shall find that the cases of real and imaginary exponents 



correspond to space-like and time-like world sheets, which both arise naturally as we shall describe. 

Af = 4 SYM is a supersymmetric gauge theory with one gauge field A M , six scalar fields 
X 1 , / = 1, 2, • ■ ■ 6 and four Weyl fermionic fields \ii & U transforming in the adjoint representation of 
the gauge group, which we take to be SU(N). The theory is conformally invariant and is specified 
by two parameters: the rank of the gauge group N and the 't Hooft coupling A, 

A = g 2 YM N. (3.1) 

(Note that the gauge coupling in the standard field theoretical convention gyM, which we shall use 
throughout, is related to that in the standard string theory convention g M by g YM = 2gj^.) 

According to the AdS/CFT correspondence, Type IIB string theory in an AdS§ x S$ spacetime 
is equivalent to an M = 4 SYM living on the boundary of the AdSs. The string coupling g s , the 
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curvature radius R of the AdS metric and the tension ^7 of the string are related to the field 
theoretic quantities as 

R 2 A 

— = v / A, 4n g a = g$ M = — . (3.2) 

a JM 

Upon first taking the large N limit at fixed A (which means g s — > 0) and then taking the large A 
limit (which means large string tension) M = 4 SYM theory is described by classical supergravity in 
AdS§ x £5. We shall describe the modification of this spacetime which corresponds to introducing 
a nonzero temperature in the gauge theory below. 

M = 4 SYM does not contain any fields in the fundamental representation of the gauge group. 
To construct the Wilson loop describing the phase associated with a particle in the fundamental 
representation, we introduce a probe D3-brane at the boundary of the AdSs and lying along n on 
S 5 , where n is a unit vector in R 6 [[45]. The D3-brane (i.e. the boundary of the AdSs) is at some 



fixed, large value of r, where r is the coordinate of the 5th dimension of AdS5, meaning that the 
space-time within the D3-brane is ordinary 3 + 1-dimensional Minkowski space. The fundamental 
"quarks" are then given by the ground states of strings originating on the boundary D3-brane and 
extending towards the center of the AdSs. 2 The corresponding Wilson loop operator has the form 

W(C) = -^TrPexp 



i (f) ds ( + n ■ XV x 2 



(3.3) 

which, in comparison with (O), also contains scalar fields X = (X , ■ ■ -X 6 ). In the large N and 



large A limits, the expectation value of a Wilson loop operator ( |3.3|) is given by the classical action 
of a string in AdS§ x S$, with the boundary condition that the string world sheet ends on the 
curve C in the probe brane. The contour C lives within the 3 + 1-dimensional Minkowski space 
defined by the D3-brane, but the string world sheet attached to it hangs "down" into the bulk of 
the curved five-dimensional AdSs spacetime. The classical string action is obtained by extremizing 
the Nambu-Goto action. More explicitly, parameterizing the two-dimensional world sheet by the 
coordinates a a = (r, a), the location of the string world sheet in the five- dimensional spacetime 
with coordinates is 

x» = x»(t,<j), (3.4) 
and the Nambu-Goto action for the string world sheet is given by 



S 

Here, 



2iia 



^—j J dadr ^/-detg af3 . (3.5) 



g a /3 = G^daX^dpx" (3.6) 



2 By the standard IR/UV connection |7!|], the boundary of the AdS^ at some large value of r corresponds to an 
ultraviolet cutoff in the field theory. The Wilson loop must be located on a D3-brane at this boundary not at some 
smaller r, in order that it describes a test quark whose size is not resolvable. Evaluating the expectation value of 
a Wilson loop then corresponds to using pointlike test quarks to probe physics in the field theory at length scales 
longer than the ultraviolet cutoff. 
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is the induced metric on the world sheet and is the metric of the 4 + 1-dimensional AdSs 
spacetime. The action ( |3.5| ) is invariant under coordinate changes of a a . This will allow us to pick 
world sheet coordinates (r, a) differently for convenience in different calculations. Upon denoting 
the action of the surface which is bounded by C and extremizes the Nambu-Goto action ( |3.5|) by 
S(C), the expectation value of the Wilson loop ( |3.3| ) is given by [15 



(W(C)) = exp[i{S(C)-S }] , (3.7) 

where the subtraction So is the action of two disjoint strings, as we shall discuss in detail below. 
To evaluate the expectation value of a Wilson loop at nonzero temperature in the gauge theory, 



one replaces AdSs by an AdS Schwarzschild black hole ||80|| . The metric of the AdS black hole 
background is given by 

r 2 1 
ds 2 = -fdt 2 + -—{dx\ + dxj + dx 2 3 ) + -dr 2 = G^dx^dx" , (3.8) 
R 2 / 

Here, r is the coordinate of the 5th dimension and the black hole horizon is at r = r^. According 
to the AdS/CFT correspondence, the temperature in the gauge theory is equal to the Hawking 
temperature in the AdS black hole, namely 

t = £. (3.io) 

The probe D3-brane at the boundary of the AdSs space lies at a fixed r which we denote r = A r . A 
can be considered a dimensionless ultraviolet cutoff in the boundary conformal field theory. We shall 
call the three spatial directions in which the D3-brane is extended x±, X2, and 23. The fundamental 
"quarks" , which are open strings ending on the probe brane, have a mass proportional to A. In order 
to correctly describe a Wilson loop in the continuum gauge theory, we must remove the ultraviolet 
cutoff by taking the A — > 00 limit. 

Now consider the set of rectangular Wilson loops shown in Fig. [l], with a short side of length L 
in the Xi-direction and a long side along a time-like direction in the t — X3 plane, which describe a 
quark- ant iquark pair moving along the x% direction with some velocity v. Here, v — corresponds 
to the loop C s tatic m Fig. [I] whereas < v < 1 corresponds to C^ t °° 1 s * cd in the figure. To analyze these 
loops, it is convenient to boost the system to the rest frame (£', x' 3 ) of the quark pair 

dt = dt' coshr/ — dx' 3 sinh^ , (3-H) 
dx% = —dt' sinh?7 + dx' 3 cosh 77 , (3-12) 

where the rapidity r\ is given by tanhr? = v. The loop is now static, but the quark-gluon plasma 
is moving with velocity v in the negative x 3 -direction. This Wilson loop can be used to describe 
the potential between two heavy quarks moving through the quark-gluon plasma or, equivalently, 
two heavy quarks at rest in a moving quark-gluon plasma "wind" . In the primed coordinates, the 
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long sides of the Wilson loop lie along t' at fixed x' 3 . We denote their lengths by T, which is the 
proper time of the quark-pair. 3 We assume that T ^> L, so that the string world sheet attached to 
the Wilson loop along the contour C can be approximated as time-translation invariant. Plugging 
( [3. 11| ) and (|3.12|) into ( |3.8|) and dropping the primes, we find 



with 



where 



r 2 1 

ds 2 = -Adt 2 - 2Bdtdx 3 + Cdx\ + — (dx\ + dx\) + -dr 2 (3.13) 

•4 = ^M). *=4|, C=£(l + 4), (3-14) 
R l \ r 4 / r l R A R l \ r 4 / 



r i = r o c °sh 2 77, T2 = r o sinh 2 77 . (3.15) 

To obtain the light-like Wilson loop along the contour Cn g ht-iike in Fig- 0> we must take the 77 — > 00 
limit. We shall see that the 77 — > 00 limit and the A — > 00 limit do not commute. And, we shall 
discover that in order to have a sensible phenomenology, we must reach the light-like Wilson loop 
by first taking the light-like limit (77 — > 00) and only then taking the Wilson loop limit (A — > 00). 
For the present, we keep both r\ and A finite. 

We parameterize the two-dimensional world sheet (|3.4j), using the coordinates 



r = t, a = x 1 E[~,^}. (3.16) 

By symmetry, we will take x M to be functions of a only and we set 

x 2 (cr) = const , x 3 (cr) = const , r = r(a) . (3-17) 
The Nambu-Goto action (|3.5|) now reads 




*=^r*J*(!*f + £). (3-18) 

with the boundary condition r(±| 1 ) = r A. This boundary condition ensures that when the string 
world sheet ends on the D3-brane located at r = vqA, it does so on the contour C which is located 
at x\ = ±4. Our task is to find r(a), the shape of the string world sheet hanging "downward in r" 



from its endpoints at r = tqA, by extremizing (p.18 ). Introducing dimensionless variables 



r = r y, a = a^, I = = nLT, (3.19) 



where T is the temperature (|3.10| ), we find that, upon dropping the tilde, 

S^VXTTf^C (3,0, 



5 In terms of the time ti a b in the rest frame of the medium, we have the standard relation T = yl — v 2 ti a b = ^j^j ■ 
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with (y 1 = d a y) 




£= t /(^_cosh 2 r / ) (1 + -^) (3.21) 



and the boundary condition y (±|) = A. In writing ( |3.20| ) we have used the fact that, by symmetry. 



y{u) is an even function. It is manifest from ( |3.20|) that all physical quantities only depend on T 
and not on R or tq separately. We must now determine y(cr) by extremizing ( |3.21| ). This can be 
thought of as a classical mechanics problem, with a the analogue of time. Since £ does not depend 
on cr explicitly, the corresponding Hamiltonian 

K^-y-f ^-^^ eonst (3.22) 

is a constant of the motion in the classical mechanics problem. 

It is worth pausing to recall how it is that the calculation of a Wilson loop in a strongly 
interacting gauge theory has been simplified to a classical mechanics problem. The large- N and 
large A limits are both crucial. Taking N — > oo at fixed A corresponds to taking the string coupling 
to zero, meaning that we can ignore the possibility of loops of string breaking off from the string 
world sheet. Then, when we furthermore take A —>■ oo, we are sending the string tension to infinity 
meaning that we can neglect fluctuations of the string world sheet. Thus, the string world sheet 
"hanging down" from the contour C takes on its classical configuration, without fluctuating or 
splitting off loops. If the contour C is a rectangle with two long sides, meaning that its ends are 
negligible compared to its middle, then finding this classical configuration is a classical mechanics 
problem no more difficult than finding the catenary curve describing a chain suspended from two 
points hanging in a gravitational field, in this case the gravitational field of the AdS Schwarzschild 
black hole. 

Let us now consider keeping A fixed and ^> 1, while increasing rj from to oo. We see that 
the quantity inside the square root in (|3.21|) changes sign when y crosses i/cosh rj. The string 



world sheet is time-like for real C (i.e. for y > ^/cosh rj) and is space-like for imaginary C (i.e. for 
y < y/coshrj). Since y = A at the boundary C, the signature of the world sheet depends on the 
relative magnitude of a/ cosh rj and A: it is time-like when yj cosh 77 < A and becomes space-like when 
\J cosh rj > A. If the world sheet in ( |3.5|) is time-like (space-like), the expectation value (|3.7|) of the 
fundamental Wilson loop is the exponent of an imaginary (real) quantity. We shall give a physical 
interpretation of this behavior in Section 3.3. Here, we explain that this behavior is consistent 
with all the phenomenology described in Section 2. For 77 = 0, the Wilson loop defines the static 
quark-antiquark potential, see (|2.2|) , and thus should and does correspond to a time-like world 
sheet. If the quark-pair is not at rest with respect to the medium, but moves with a small velocity 
v = tanh rj, one still expects that the quark-pair remains bound and the world-sheet action remains 
time-like. We shall see, however, that for large enough 77 a bound quark-antiquark state cannot 
exist. Once we reach rj = 00, namely the light-like Wilson loop which we saw in Section 2 originates 
from eikonal propagation in high energy scattering and is relevant to deep inelastic scattering, the 
Cronin effect, and jet quenching, in order to have a sensible description of these phenomena we see 
from ( 2.15|) or equivalently (|2.6|) that the expectation value of the Wilson loop must be the exponent 



- 16 - 



of a real quantity. This expectation is met by ( |3.21| ) since the string world sheet is space-like as 
long as y/ cosh 77 > A. This demonstrates that in order to sensibly describe any of the applications 
of Wilson loops to high energy propagation, including in particular in our nonzero temperature 
context the calculation of the jet quenching parameter q, we must take 77 — > 00 first, before taking 
the A — > 00 limit. 

In subsection 3.1, we shall review the calculation of the quark-antiquark potential and screening 
length as a function of the velocity v. In subsection 3.2, we calculate the jet quenching parameter. 
And, in subsection 3.3, we return to the distinction between the time-like string world sheet of sub- 
section 3.1 and the space-like string world sheet of subsection 3.2, and give a physical interpretation 
of this discontinuity. 

3.1 Velocity-dependent quark-antiquark potential and screening length 

In this subsection we compute the expectation values of Wilson loops for \J cosh 17 < A, from which 
we extract the velocity-dependent quark-antiquark potential and screening length. At the end of 
the calculation we take the heavy quark limit A — ► 00. In fact, because we are interested in the 
case a/ cosh rj < A, in this subsection we could safely take A — > 00 from the beginning. The results 
reviewed in this subsection were obtained in Refs. |3"5| . 

We denote the constant of the motion identified in Eq, (|3.22|) by q, and rewrite this equation 

as 

y' = - q ^(y 4 -i)(y A -yt) (3.23) 

with 

y\ = cosh 2 77 + g 2 . (3.24) 

Note that y 4 > cosh 2 77 > 1. The extremal string world sheet begins at a = — 1/2 where y = A, and 
"descends" in y until it reaches a turning point, namely the largest value of y at which y' = 0. It 
then "ascends" from the turning point to its end point at o = +£/2 where y = A. By symmetry, the 
turning point must occur at o = 0. We see from (|3.23| ) that in this case, the turning point occurs 



at y = y c meaning that the extremal surface stretches between y c and A. The integration constant 
q can then be determined 4 from the equation ~ = J Q 2 da which, upon using ( |3.23| ), becomes 



l = 2q dy ; (3.25) 

4 VG/ 4 - - 1) 

The action for the extremal surface can be found by substituting ( |3.23| ) into ( |3.20| ) and ( |3.21| ), 

S(t) = rxTrf\ / t ~^ (3.26) 



Equation ( |3.26|) contains not only the potential between the quark-antiquark pair but also the 



static mass of the quark and antiquark considered separately in the moving medium. (Recall that 



'For equation ( |3.23| ) to be well denned, we need < q 4 < A 4 — cosh 2 rj. 
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we have boosted to the rest frame of the quark and antiquark, meaning that the quark-gluon plasma 
is moving.) Since we are only interested in the quark- antiquark potential, we need to subtract from 
( p.26| ) the action Sq of two independent quarks, namely 



E{L)T = S(l) — S 



(3.27) 



where E(L) is the quark-antiquark potential in the dipole rest frame. The string configuration 
corresponding to a single quark at rest in a moving hot medium in M = 4 SYM was found in 
Refs. 123, |24l , from which one finds that 



S = VXTT J dy. 



(3.28) 



To be self-contained, in appendix |A] we review the solution of [23, 24], along with a family of new 
drag solutions describing string configurations corresponding to mesons made from a heavy and a 
light quark. 

To extract the quark-antiquark potential, we use ( P-25| ) to solve for q in terms of I and then 
plug the corresponding q{l) into ( |3.26| ) and (|3.27|) to obtain E(L). We can safely take the A — > oo 
limit, and do so in all results we present. We show results at a selection of velocities in Fig. In 
the remainder of this subsection, we describe general features of these results. 

First, Eq. (|3.26|) has no solution when / > l ma ^(r]), where Z max is the maximum of l(q). We see 
that / max decreases with increasing velocity. 

We see from the left panel in Fig. that for a given I < Z ma x(^)> there are two branches of 
solutions. The branch with the bigger value of q, and therefore the larger turning point y c , has 
the smaller E(L) — corresponding to the lower branches of each of the curves in the right panel 
of the figure. The upper branches of each curve correspond to the solutions for a given I < Z max 
with smaller q and y c . Because they have higher energy, it is natural to expect that they describe 
unstable solutions sitting at a saddle point in configuration space fl32| , [34|[ . This has been confirmed 
explicitly in Ref. 
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When T] is greater than some critical value r] c , E(L) is negative for the whole upper branch. 
When r) < r)c, there exists a value l c (rj) < l max such that the upper branch has an E(L) which is 
negative for I < l c and positive for I > l c . l c goes to zero as r\ goes to zero. If r\ < r\ c and I > l c , 
then if the unstable upper branch configuration is perturbed, after some time it could settle down 
either to the lower branch solution or to two isolated strings each described by the drag solution 



of Ref. [^, |24[ and Appendix A. (Note that E > means that a configuration has more energy 
than two isolated strings.) On the other hand, if E(L) is negative for the upper branch, when this 
unstable configuration is perturbed, the only static solution we know of to which it can settle after 
some time is the lower branch solution. 

We see from Fig. |2| that using the action of the dragging string solution of Refs. |23, |24| as So as 
we do and as was considered as an option in Ref. |3~5(|, ensures that the small-distance behavior of the 



potential is velocity- independent. This seems to us a physically reasonable subtraction condition; 
it is analogous to the renormalization criterion used to define the quark-antiquark potential in 
lattice calculations, namely that at short distances it must be medium-independent f43fl. Choosing 
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0.4 0.6 

q 1 = 7ULT 

Figure 2: Left panel: the quark-antiquark separation l(q) as a function of the integration constant q for 
a quark-antiquark dipole oriented orthogonal to the wind propagating at different velocities v = tanhr/. 
We discuss the case where the dipole is not orthogonal to the wind in Section 5. Right panel: The qq 
static potential for the same quark-antiquark configurations as in the left panel. Note that the potential 
is normalized such that the small-distance behavior of the potential is unaffected by velocity-dependent 
medium effects. 



the velocity-dependent subtraction (|A.12|) instead, considered as an option in Ref. p5| , makes the 
unstable upper branch have lirn^o E(L) = for all velocities, but in so doing makes the stable 
lower-branch have a velocity-dependent E(L) at all L, including small L. 



One can obtain an analytical expression for l r 
in powers of gives 



in the limit of high velocity. Expanding ([3.25 ) 



l(q) 




(3.29) 
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Truncating this expression after the second term, we find for the maximum 

I = ^ T ® (_?_ + 1 + .. 

max 3«/*r(J) Vcosh 1 / 2 ^ 5 cosh 5 / V 

= 0.74333 ( + — I J . (3.30) 

Vcosh 1/2 r/ 10cosh 5/2 r? / 

Note that L max = %p can be interpreted as the screening length in the medium, beyond which the 
only solution is the trivial solution corresponding to two disjoint world sheets and thus E(L) = 0. 
The first term of this expression was given in |^TJ (see also ||32||), the second term in [j35|l . As we 
shall discuss further at the end of Section 5, if we set rj = in ( |3.30| ), this expression which was 
derived for rj —>■ oo is not too far off the rj = result, which is £ max = 0.869. Hence, as discovered in 
Ref. [31 1, the screening length decreases with increasing velocity to a good approximation according 



to the scaling 

L max (v) ~ = (3.31) 

cosh 1 7] V7 



with 7 = 1/ yl — v 2 . This velocity dependence suggests that L s should be thought of as, to a good 
approximation, proportional to (energy density) _1 / 4 ; since the energy density increases like 7 2 as 
the wind velocity is boosted. 

If the velocity-scaling of L s that we have discovered holds for QCD, it will have qualitative con- 
sequences for quarkonium suppression in heavy ion collisions |31"|| . For illustrative purposes, consider 



the explanation of the J/^f suppression seen at SPS and RHIC energies proposed in Refs. [^T], j82| : 
lattice calculations of the gg-potential indicate that the J/^(1S) state dissociates at a temperature 
~ 2.1T C whereas the excited Xc(2P) and ^'(2S) states cannot survive above ~ 1.2T C ; so, if collisions 
at both the SPS and RHIC reach temperatures above 1.2T C but not above 2.1T C , the experimental 
facts (comparable anomalous suppression of J/^f production at the SPS and RHIC) can be under- 
stood as the complete loss of the "secondary" J/^'s that would have arisen from the decays of the 
excited states, with no suppression at all of J/^'s that originate as J/^'s. Taking Eq. ( |3.31| ) at 
face value, the temperature T diss needed to dissociate the J/ty decreases oc (1 — f 2 ) 1 / 4 . This indi- 
cates that J/ty suppression at RHIC may increase markedly (as the J /iff (IS) mesons themselves 
dissociate) for J/^'s with transverse momentum pt above some threshold that is at most ~ 9 GeV 
and would be ~ 5 GeV if the temperatures reached at RHIC are ~ 1.5T C . The kinematical range in 
which this novel quarkonium suppression mechanism is operational lies within experimental reach 
of future high-luminosity runs at RHIC and will be studied thoroughly at the LHC in both the J/^f 
and Upsilon channels. If the temperature of the medium produced in LHC collisions proves to be 
large enough that the J/^(1S) mesons dissociate already at low p?, the p^-dependent pattern that 
the velocity scaling ( |3.31| ) predicts in the J/\I/ channel at RHIC should be visible in the Upsilon 
channel at the LHC. 

As a caveat, we add that in modelling quarkonium production and suppression versus pt in 
heavy ion collisions, various other effects remain to be quantified. For instance, secondary pro- 
duction mechanisms such as recombination may contribute significantly to the J/^> yield at low 
Pt, although the understanding of such contributions is currently mo del- dependent. Also, at very 
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high Pt, J/y& mesons could form outside the hot medium [83|. Parametric estimates of this effect 



suggest that it is important only at much higher pt than is of interest to us, and we are not aware 
of model studies which have been done that would allow one to go beyond parametric estimates. 
The quantitative importance of these and other effects may vary significantly, depending on details 
of their model implementation. In contrast, Eq. ( |3.31| ) was obtained directly from a field-theoretic 
calculation and its implementation will not introduce additional model-dependent uncertainties. 
For this reason, the velocity scaling established here must be included in all future model calcula- 
tions. We expect that its effect is most prominent at intermediate transverse momentum, where 
contributions from secondary production die out or can be controlled, while the formation time 
of the heavy bound states is still short enough to ensure that they would be produced within the 
medium if the screening by the medium permits. 

3.2 Light-like Wilson loop and the jet quenching parameter 

In order to calculate the jet quenching parameter we need to take the rj — > oo limit in which the 
Wilson loop becomes light-like first, with the location of the boundary D3-brane A large and fixed, 
and only later take A — > oo. As we approach the light-like limit, it is necessary that ^cosh r\ > A. 
In this regime, as we discussed below equation ( |3.22| ), the world sheet is space-like, meaning that 
the expectation value of the Wilson loop is the exponential of a real quantity. As we reviewed in 
Section 2, this must be the case in order to obtain sensible results for both medium-induced gluon 
radiation of Eq. ( |2.8|) and the virtual photo-absorption cross section in deep inleastic scattering of 
Eq. O- 

When y/ cosh 77 > A, the first order equation of motion, given by (|3.22|) , reads 



1 

Q 

with 



y ,2 = - 2 {y A -^){ytn-y') (3.32) 



Vm = cosh 2 - q 2 . (3.33) 

The consistency of ( |3.32| ) requires that y m > A, which implies that the integration constant q is 
constrained to < q 2 < cosh 2 rj — A 4 . Equation (|3.32|) has a trivial solution 

y(p) = A = const, q 2 = cosh 2 77 - A 4 . (3.34) 



However, one can check that ( |3.34|) does not solve the second order Euler-Lagrange equation of 



motion derived from ( |3.21p and thus should be discarded. Because y m > A, the nontrivial solution 
of ( |3.32|) which descends from y = A at a = —1/2 descends all the way to y — 1, where y' = 0. 
Thus, for any value of I the string starts at y = A and descends all the way to the horizon, where 

it turns around and then ascends back up to y = A. The integration constant q can be determined 

1 

from the equation ^ = J Q 2 da, i.e 

,A x 

l = 2q dy - (3.35) 
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upon using ( |3.32|) . The action (|3.18| ) takes the form 

S(l) = i^XTT f dy 7 h ' " 7 4 . (3.36) 

This action is imaginary and corresponds to a space-like world sheet. 

To extract q introduced in ( |2.6| ), ( [2.11| ) and ( |2.15| ), we first take rj — > oo, making the contour C 
light-like, and only then take the A — > oo limit needed to ensure that we are evaluating W(C), with 
the end of the string on the D3-brane at y = A following the contour C precisely, q can be obtained 
by studying the small /-dependence of the action (|3.36|) , which can be done analytically. We start 
from the expansion of (|3.35|) , 

I = T dy^=L= + o(q 3 , -^r-) . (3.37) 

cosh rj J x ^Jy± - 1 V cosher// 



Upon defining 



A^oo 



A 



: r(!) 



a = lim / dy — -j^=^= = y/n — )4hr , (3.38) 



Vv^i r(i)' 



we find that in the small / (equivalently, small q) limit 

I = ^- . (3.39) 
cosh 7/ 

In the same limit, the action (|3.36|) takes the form 

S(l) = S ,(0) + g 2 S (1) + O(g 4 ), (3.40) 



where 



S^= l VXTT dy.rV 1 ^ 1 , (3-41) 



q 2 S^(l) 




(cosh 2 r] - y 4 ) {y 4 -l) 
iy/XTTq 2 a .V\tt 2 T 3 



2 cosh rj 8 a 



(T cosh 7]) L , (3.42) 



where we have used ( p. 38 ), ( |3.39|) and / = it LT . Also, we have kept the dominant large r]- 
dependence only. We identify (Tcoshri) = L~ / \/2, where L~ is the extension of the Wilson loop 
in the light-like direction, entering in ( [2.11|) and ( |2.15|) . 

As in Section 3.1, in order to determine the expectation of the Wilson line we need to subtract 
the action of two independent single quarks, this time moving at the speed of light. In Appendix [A], 
we analyze the string configuration corresponding a single quark moving at the speed of light. There 
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we find a class of solutions with space-like world sheets and also a class of solutions with time-like 
world sheet. Our criterion to determine which solution to subtract is motivated from the physical 
expectation discussed in Section 2, i.e. 

lim \S(l) - S ] = 5 (0) -£0 = 0. (3.43) 

z->o 

Among the classes of solutions discussed in Appendix A, the only one satisfying ( |3.43| ) is the space- 
like world sheet described by Eqs. (|A.17|) and ( |A.18| ) with p — 0. In this configuration, So is the 
action of two straight strings extending from y = A to y = 1 along the radial direction and is given 
by 

So = i VXTTj\yJ^^. (3.44) 



The L 2 -term in the exponent of ( |2.15| ) can then be identified with the 0(L 2 )-term ( |3.42 ) of the 



action S(l), and we thus conclude that the jet quenching parameter in (|2.15|) is given by 

vr^r (§) 



qSYM 



r(f) 



-V\T 3 . (3.45) 



We have used the fact that, as in ( |2.9| ), in the large- N limit the expectation value of the adjoint 
Wilson loop which defines q in (|2.15|) differs from that of the fundamental Wilson loop which we 



have calculated by a factor of 2 in the exponent S. 

In Ref. [)15| , the result (|3.45|) was obtained starting directly from the loop Cught-iike, described 



in the rest frame of the medium using light-cone coordinates. Here, we showed that one can obtain 
the same result by taking the v — > 1 limit of a time-like Wilson loop. It is also easy to check that 
the trivial solution ( |3.34j ) goes over to the constant solution discussed in ||15|| , which has a smaller 



action than ( |3.36|) . In [15| this trivial solution was discarded on physical grounds. Here, we see 



that if we treat the light-like Wilson line as the r] — > 00 limit of a time-like one, this trivial solution 
does not even arise. We also note that in the light-like limit, the coefficient in front of the scalar 
field term in ( [3. 3D goes to zero and (|3.3|) coincides with (|2.1|) . 

In Section 4 we shall determine how q changes if the medium in which the expectation value of 
the light-like Wilson loop is evaluated has some flow velocity at an arbitrary angle with respect to 
the direction of the Wilson loop. In Section 6 we shall discuss the comparison between our result 
for q and that extracted by comparison with RHIC data, as well as discuss how q changes with the 
number of degrees of freedom in the theory. 

3.3 Discussion: time-like versus space-like world sheets 

We have seen that as we increase 77 from to 00 while keeping A fixed and large, the behavior of 
the string world sheet has a discontinuity at ^ s /coshr] = A, below (above) which the world sheet is 
time-like (space-like). Here we give a physical interpretation for this discontinuity. Recall first from 
Section 3.1 that if cosh 77 ^> 1 but \J cosh r] < A, the screening length L max is given by 

^max = ^L- . (3.46) 

7TV cosh r\ 1 
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Next, note that the size 5 of our external quark on the D3-brane at y = A, i.e. at r = r A = tcR 2 TA 
can be estimated using the standard IR/UV connection, namely ]75| 



5 ~ ^ ~ , (3-47) 
M AT v ; 

where M = ±\/ATA is the mass of an external quark as can be read from (|3.28|) . (The apparent 
T-dependence of Q3.47| ) is due to our definition of A, with the ultraviolet cutoff given by Ar , and 
does not reflect genuine temperature-dependence.) Thus the condition y/coshrj = A corresponds 



to 32 



S . (3.48) 



When a/ cosh rj <C A, meaning that 5 <C L max , we expect that if instead of merely analyzing Wilson 
loops we were to actually study mesons, we would in fact find a bound state of a quark and anti- 
quark. In this regime, it is reasonable to expect that the expectation value of the Wilson loop should 
yield information about the quark-antiquark potential, meaning that it must be the exponential of 
an imaginary quantity meaning that the string world sheet must be time-like, as indeed we find. 
On the other hand, when y/coshrj 3> A, meaning that 5 ^> L max , the size of one quark by itself is 
much greater than the putative screening length. This means that the quark and antiquark cannot 
bind for any L, meaning that the transition at y/ cosh 77 ~ can be thought of 

"deconfinement" or "dissociation" transition for quarkonium mesons made from quarks with mass 
M. Furthermore, in a regime in which the size of one quark is greater than the putative screening 
length, the concept of a quark-antiquark potential (and a screening length) makes no sense. Instead, 
in this regime it is appropriate to think of the quark-antiquark pair as a component of the wave 
function of a virtual photon in deep inelastic scattering, and hence to think of the Wilson loop as 
arising in the eikonal approximation to this high energy scattering process, as discussed in Section 
2. From our discussion there, it is then natural to expect a space-like world sheet, which gives the 
desired (W) ~ exp[— S] behavior with S real. 

Our discussion explains the qualitative change in physics, but it does not explain the sharpness 
of the discontinuity that we find at yj cosh 77 = A, which likely has to do with the classical string 
approximation (which corresponds to large N and large A limit) we are using. When a/ cosh 77 < A 
there is a discontinuity between L < L max and L > L max where the quark-antiquark potential goes 
from being nonzero to zero. This discontinuity is smoothed out by finite A corrections, with the 
exponentially small quark-antiquark potential at large distances corresponding to physics that is 
nonperturbative in a'. Presumably the discontinuity at y/ cosh rj = A is also smoothed out at finite 
A and N. Further insight into this question could perhaps be obtained without relaxing the large- N 
and large-A limits by studying mesons rather than Wilson loops. 



The operational consequences of the discontinuity at \J cosh 77 = A are clear. To compute the 
quark-antiquark potential and the screening length in a moving medium, we take A to infinity at 
fixed 77. To compute q, we must instead first take the 77 — > 00 limit at finite A, and only then take 
A — > 00. The two limits do not commute. 
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4. The jet quenching parameter in a flowing medium 

In section [T2], we have evaluated the expectation value of a light-like Wilson loop specified by the 
trajectory of a dipole moving in a light-like direction, k 11 = (l,n) with n a unit vector and hence 
k 2 = 0. The world sheet defined by this light-like Wilson loop is space-like and the behavior of 
the Nambu-Goto action in the limit of small dipole size determines the jet quenching parameter 
q. It has been argued previously that the motion of the medium orthogonal to the trajectory of 
the dipole can affect the value of q in a nontrivial fashion [85]. Furthermore, if the medium is 
flowing parallel to or antiparallel to the trajectory of the dipole with velocity Vf = tanh rjf, there is a 
straightforward effect on q: the calculation goes through unchanged, with L~ understood to be the 
light-cone distance in the rest frame of the medium, but the relation between L~ and the distance 
Az travelled in the lab frame is modified: Az = (L~/y/2) exp(^j), where the sign convention is 
such that rjf > corresponds to the dipole velocity and flow velocity parallel (i.e. the dipole feels a 
"tail wind") while rjf < means that the dipole feels a head wind. Correspondingly, q is multiplied 
by a factor of exp(— rjf), meaning that it increases in a head wind and decreases in a tail wind. In 
this Section, we calculate how the jet quenching parameter q depends on the speed and direction 
of the collective flow of the medium, allowing for any angle between the jet direction and the flow 
direction. 

The calculation of the effect on jet quenching parameter q due to the collective motion of the 
medium turns out to be straightforward, once the geometry of the problem is set up. We shall 
specify the light-like four-momentum k^ = (l,n) (the direction of motion of the hard parton which 
is losing energy; the direction of propagation of the dipole moving at the speed of light which defines 
the Wilson loop) by taking n to point along the negative ^-direction. According to the way the 
BDMPS energy loss calculation is set up, the dipole is always perpendicular to the direction of its 
motion, so we choose the dipole orientation to point in a direction rh which must lie in the (x\, 22)- 
plane. Now, we set the medium in motion. The most general "wind velocity" has components 
parallel to and orthogonal to the dipole direction n. We choose v — vl to lie in the (X2, X3)-plane. 
Because we fix the orthogonal component to lie along the 22-direction, we must leave the direction 
of the dipole orientation in in the (xi, a^-plane unspecified. Thus the most general configuration is 
described by four parameters, the transverse separation L of the Wilson loop in the lab frame and 

coshr] f = 1 9 = Z(l,h), <p = Z(rh,x 1 ). (4.1) 
V 1 — v 2 

In the lab frame, the trajectory of the end points of the dipole can be written as 

A± = m ± , (4.2) 

where 

fc" = (l,7i), m" = (0,m), k-m = 0, k 2 = 0, m 2 = 1 . (4.3) 
Now, we boost A 11 with v = vl, boosting into a frame in which the medium is at rest. We obtain 

A± = k'H ± -m'v , (4.4) 
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for some k'^ and m' M which again satisfy 



k = 0, m' 2 = 1, rri ■ k' = . (4.5) 

In general, m ,fl has a nonzero 0-th component and thus the two ends of the dipole do not have the 
same time. To fix this we write 

Af = k'»(t±t )±^m"», (4.6) 

where we have defined 

m"" = m'» - tok' 11 , (4.7) 

and choose t such that the zeroth component of m" is zero, making m"^ purely spatial. It is easy 
to confirm that, given ( |4.6| ), we now have 

k /2 = 0, m ,,2 = l, m ".k' = 0. (4.8) 

We now have almost exactly the same Wilson loop configuration as we had in our original calculation 
of Section 3.2 when the medium was at rest from the beginning, with the only difference being that 
the two long sides of the Wilson loop do not start and end at equal times, due to the shift £ - 
This is immaterial when L_ is big: in our evaluation of the Wilson loop we always assumed time 
translational invariance anyway, neglecting the contribution of the "ends of the loop" relative to 
that of the long, time translation invariant, mid-section of the loop. We thus find that in the 
presence of a wind velocity 

(W A (C)) = exp [S(C)} (4.9) 

with 

S(C) = -J-h(ir)'L 2 , (4.10) 

where go is the value with no wind and where 

(L-)' = V2k'°t = k'°L- (4.11) 

is the light-cone distance travelled in the rest frame of the medium whereas L~ is the corresponding 
quantity in the lab frame. We thus conclude that the only effect of the collective flow of the 
medium on q is what we called the straightforward effect above, namely that due to the Lorentz 
transformation of L~ . From the standard Lorentz transformation rule, 

k'° = cosh i] f — sinh T]f(l • n) = cosh rjf — sinh rjf cos 9 . (4-12) 

We thus find 

q = (cosh rjf — sinh rjf cos 6) q . (4-13) 

This result is independent of 0. 

We have established the transformation rule ( 4 . 1 3|) by boosting to the rest frame of the medium. 



This reduced the problem to one with no wind but with a Lorentz transformed longitudinal ex- 
tension ( [LI If ). Alternatively, the same result ( [4.13 ) can be obtained by starting from the metric 
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corresponding to the medium having a velocity v and doing the Wilson loop computation in this 
metric. We have confirmed by explicit calculation for several examples that the same result ( OB ) 
is obtained. 

The derivation of the scaling ( [4.13[ ) relied only on properties of Lorentz transformations; nothing 
in the calculation of the underlying q (which depends on T and iV and A in M = 4 SYM and varies 
from one theory to the next as we shall discuss in Section 6) comes in. We conclude that the scaling 
( [4.13| ), which describes how the jet quenching parameter q depends on the collective flow velocity 
of the medium doing the quenching, applies in QCD also. R. Baier et al. have reached the same 
conclusion independently [ |S6|j . 

To get a sense of the order of magnitude of the effect, we note that transverse flow velocities 
in excess of half the speed of light are generated by the time the matter produced in a heavy ion 
collision freezes out. A velocity v = 0.5 corresponds to rjf = 0.549, which yields q = 1.732 q for a 
head wind (8 = tt), q = 1.155 go for 9 = n/2, and q = 0.577 for a tail wind (6 = 0). An investigation 
of the quantitative consequences of (|4.13|) requires modelling of the geometry and time-development 
of the collective flow in a heavy ion collision, along the lines of the analysis in Refs. 



5. The static qq potential for all dipole orientations with respect to the 
wind 



In section |3.1] , we analyzed the quark-antiquark potential for a gg-dipole which was oriented in the 
Xi direction and which propagated orthogonal to its orientation along the x% direction with velocity 
v = tanh^. Here, we extend this analysis to the case where the dipole is tilted by an arbitrary angle 
9 with respect to its direction of motion, see Fig. 13. For 6 = ir/2, we recover the results obtained 



in section 3.1 above 



We work in the boosted metric ( |3.13| ), in which the dipole is at rest. The dipole lies in the 



(xi, X3)-plane and the parametrization of the two-dimensional world sheet of the corresponding 
Wilson loop is 

r = t , a = X\ , X2 = const. , X3 = £3(0") , r = r(cr) . (5.1) 

For a dipole with length L whose orientation makes an angle 9 with its direction of propagation (the 
^-direction) the projections of the dipole on the X\ and X3 axis are of length Lsin^ and Lcos8, 
respectively. We define dimensionless coordinates 

V = ~, z = x 3 -, a = a-, l = L-, (5.2) 
and drop the tilde. The boundary conditions on y(a) and z(a) then become 

y ( ±- sin 6 ) = A , z ( ±- sin 6 ) = ±- cos 6 . (5.3) 



2 \ 2 2 
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Figure 3: Schematic picture of a gg-dipole, moving along the X3-direction. The dipole is oriented along 
an arbitrary orientation 9 in the (x\, X3)-plane. The trajectories, along which the quark and antiquark 
propagate, specify the boundary C of a two-dimensional world sheet, which extends into the 5-th bulk 
dimension y = r/r$. The shape of this world sheet is characterized by the functions y(cr) and z(a), where 
a = x\. By symmetry, y has a turning point at y(a = 0) = y c . For generic values of the angle 6, the 
solution z(a) deviates from a straight line. 



Following the calculation of section |3.1| , the Nambu-Goto action for (|5.1| ) can be written in the form 

ril 
hi 



( ggTj ), namely S(C) = V\TT j l J 2 daC, with the La grangian reading 



C=J (y* - cosh 2 V ) 1 1 + ) + z>* (y4 - 1) , (5.4) 



where y' and z' denote derivatives with respect to a . The Hamiltonian is 



H - c - v'— - z'— - y'- cosh ^ - a (55) 
H ~ L V dy> Z dz>~ C ~ q > 

a constant of the motion. The momentum conjugate to z 

^> = — z=v (5 - 6) 

is also a constant of the motion. For a time-like world sheet, the constants of motion q and p must 
be real. The equations of motion can be written in the form 

q 2 y' 2 = (y 4 - cosh 2 ^ 4 - 1 -p 2 ) - q 2 {y' - 1) , (5.7) 

2/2 2 ( V 4 -cosh 2 rA 2 

q z = p y4 _ 1 J ■ (5.8) 
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Figure 4: Solutions y(cr), z{p) of the differential equations (5.7) and (|5.8|) with boundary conditions (|5.3[) 
for different orientation angles and rapidity rj. These solutions characterize the two-dimensional world 
sheet bounded by the Wilson loop. For definitions, see Fig. || and text. Left panel: y(a) for rj = 1 for 
varying 9. Middle panel: z(a) — a /tan 9 for rj = 1 for varying 9. As Fig. || illustrates, z(a) — a /tan 6 
would be zero if the projection of the string onto the D3-brane at y = A were a straight line. Third panel: 
z(a) — aj tan 9 for 9 = 7r/4 for varying 77. 



Generic features of their solutions have been pointed out in Ref. [aT| already. Fig. [| shows numerical 



results. Since y' depends only on y and since the boundary condition (|5.3|) for y is symmetric under 
o — > — cr, y(cr) must be an even function of a. It descends (y' < 0) for — l/2sm6 < o < and then 
ascends for < a < 1/2 sin8. These features are clearly seen in Fig. f|. The turning point y c = y(0) 
satisfies the condition 

(yt - cosh 2 77) (yt - 1 - p 2 ) - g 2 {y 4 c - l) = . (5.9) 

Connecting the gg-pair by a straight line in the (xi, X3)-plane would correspond to 2(0") = a/ tan#. 
To test for deviations of the string world sheet away from this straight line, we plot z{a) — a j tan# 
in Fig. ^. We find a deviation of sinusoidal form for all angles except 9 = 0,7r/2. As an aside, we 
note that if one thinks of the two-dimensional world sheet as a flat piece of paper, draws on it a 



-29- 



o 
o 
o 



o 
o 



& 2 



3 



— e = 7t/i80 

— e = 7i/6 
e = 7t/3 



0.01 



100 



Figure 5: The integration constant p as a function of the integration constant q for rj = 1 and several 
fixed values of the orientation angle 9. This relation is defined by eqs. (5.10) and (5.11). Along one of 
these curves, I changes with q as shown in Fig. |(| 



straight line connecting q and q, and rolls it up as depicted in Fig. |3], then the projection of this 
straight line on the (x%, a^-plane would show a qualitatively similar sinusoidal wiggle. However, 
the use of this analogy is limited, since we cannot specify in which sense or to what extent the 
two-dimensional world sheet is flat. Also, the observed deviation from the straight line behavior 
z(a) = o '/ tan 9 depends on rapidity. For rj = 0, no deviation is possible since no direction in the 
(xi, x 3 )-plane is singled out. For increasing values of 77, the deviation increases as seen in Fig. [|. 

The constants q and p must be related to the values of I and 6. The relationships are obtained 
by integrating the equations of motion (|5.7|) and (|5.8|), giving 



- sin 9 
2 



dy 



(y 4 - cosh 2 rj) (y 4 - 1 - p 2 ) - q 2 (y 4 - 1) 



(5.10) 



- cos 9 = p 
2 y 



y 4 — cosh r\ 



dy 



y 4 -l 



(y 4 - cosh 2 rj) {y 4 - 1 - p 2 ) - q 2 {y 4 - 1) 



(5.11) 
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With q and p determined, the qq static potential in a moving thermal background then reads 



E(l,9; V )=S(C)-S 



VXTT 



lie 



(y A — cosh 2 7]) 



(y 4 - cosh 2 rj) (y 4 - 1 - p 2 ) - q 2 (y 4 - 1) 



- 1 



dy 



V\TT (y c 



(5.12) 



Here, the subtraction term So, given in (|3.28|) is the action for two isolated strings described by the 



dragging solution of Refs. [23, 23 1 and Appendix A. 

We have evaluated the potential E(l, 9; rj) as a function of the size I of the dipole, its orientation 
9 with respect to its direction of motion, and its velocity v = tanhr/ with respect to the thermal 
heat bath. Since the potential in ( |5.12[) is written in terms of the integration constants q and p, 
it is useful to determine first how p depends on q for fixed 9 and rj. To do this, we write tan# as 
the ratio of Eqs. (|5.10|) and ( 5.11 ). We find p(q) to be a monotonously increasing function, whose 
slope decreases with increasing angle 9, see Fig. |5|. For the maximal angle 9 = tt/2, p vanishes 
independent of the value of q. This is the case of a dipole oriented orthogonal to the wind, where 
Eq. (|5.10| ) reduces to Eq. ( p. 251 ), an d the present calculation becomes that of Section [3~T1 . For the 
opposite limit of a dipole oriented parallel to the wind, 9 = 0, the parametrization (5TT), (|5T 



of the 

two-dimensional world sheet does not apply. However, the Nambu-Goto action is reparametrization 
invariant and, as described in Appendix |B|, in a parametrization which is suitable for < 9 < n/2 
we find that the measurable quantity E(l, 9, rf) depends smoothly on 9 for 9 — > 0. 

Knowing p(q) for fixed r] and 9, the rescaled dipole size l(q) = l(q,p(q)) can be written as a 
function of q only. It takes values in the range I G [0, / max ]. Here, the maximal dipole size / max is 
the screening length above which bound states do not exist. The value of q at which the maximum 
of l(q) occurs depends strongly on the angle 9, as shown in the left panel of Fig. This is a feature 
of our parametrization: for smaller angle 9, p(q) is a more steeply rising function (see Fig. ||), and 
most of the ^-dependence of l(q,p(q)) comes from the p(g)-dependence. The value of Z max decreases 
slightly with increasing angle 9. This is consistent with the expectation that the dipole is easier to 
dissociate if it is oriented orthogonal to the direction of the wind, but the effect is slight. 

With p(q) determined, the gg-static potential E(l(q,p(q)),9,rj) also becomes a function of q 
only. This defines curves {l(q), E(l(q))}, parametrized by the integration constant q G [0, 00]. The 
qq static potential (|5.12|) is a double- valued function of I in the range / G [0, Z max ], see the right panel 
of Fig. ^. The configurations whose energy is given by the upper branch of E(l) are presumably 
unstable, as has been shown explicitly for 9 = n/2 in Ref. [p6 |. The lower branch displays the 
typical short-distance behavior of a qq binding potential. For fixed rapidity r\, this potential shows 
the expected ^-dependence: the qq pair is more strongly bound if the dipole is aligned with the 
direction of motion, and this binding decreases as the dipole presents itself at a larger angle with 
respect to the wind, see Fig. ||. 

In Fig. ^| in Section 3, we have explored the ^-dependence of the qq static potential for a dipole 
oriented orthogonal to the wind. The screening length displays the dominant Lorentz-7 dependence 



I miB oc 1/ y/ cosh r] = 1/y/j, as given in ( |3.30|) . This velocity dependence is much stronger than the 
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Figure 6: Left panel: The size I of the dipole and its orientation angle 9 define the integration constants 
q and p, see eqs. fTToD and ( fTllD . The plot shows I as a function of q for several fixed angles 6. Along 
each of these curves p varies as shown in Fig. |5[ Right panel: The static g ^-potential ( |5.12 ) for rapidity 
T) = 1 and different orientation angles of the dipole with respect to the direction of motion. 



angular dependence displayed in Fig. || The velocity dependence of the potential E(l), shown on the 
right hand side of Fig. [], also shows clearly that the short distance behavior of the potential is not 
affected by velocity-dependent medium effects. This is a consequence of choosing the regularization 
prescription ( |3.28| ). 

Finally, we show in Fig. |7| the screening length l max multiplied by y^coshrj. We include curves 
for 9 = and 9 = 7r/2; those for other angles lie in between these two. Note that both curves have 
the same value of Z max in the rj — > limit as they must. The flat behavior of these curves at large r\ 
illustrates that Z max oc 1/ \J cosh r\ is the leading large-?] dependence for all dipole orientations. This 
leading behavior provides a numerically very accurate approximation (< 1% deviation) for 77 > 2, 
and even for 77 = 0, it is accurate to within 20% (note the suppressed zero in Fig. 0). Including the 
(coshr?) -5 / 2 term in the analytical expansion (|3.30| ) improves the description. 
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Figure 7: The screening length / max times its leading large-77 dependence U cosh(r/) . The exact results 
are given for dipoles oriented perpendicular to the wind {6 = vr/2) and parallel to the wind (0 = 0). The 
9 = it/2 curve is compared to the analytical large-?] approximation (|3.30| ). Keeping only the first term in 
this analytical expression corresponds to a horizontal line on the figure; including the term proportional 
to (coshr?)" 5 / 2 improves the agreement with the exact result. 



6. Discussions and Conclusions 



In Section 2, we reviewed the physical arguments why the expectation value of the time-like Wilson 
loop which describes the quark-antiquark potential must be the exponential of an imaginary quantity 
whereas the expectation value of the light-like Wilson loop which arises in the physics of deep 
inelastic scattering, proton-nucleus collisions, and the calculation of the jet quenching parameter 
relevant to parton energy loss in heavy ion collisions is instead the exponential of a real quantity. 
In Section 3, we saw how these results emerge by direct calculation in M = 4 SYM theory at strong 
coupling, where via the AdS/CFT correspondence the calculation of the expectation values of these 
two types of Wilson loops reduces to the evaluation of the action of an extremal string world sheet, 
time-like in the first case and space-like in the second. These aspects of our paper are discussed at 
length in Section 3 and we shall not discuss them further here. 

This section incorporates several different discussions, while along the way summarizing many 
of our conclusions. In Section 6.1, we summarize what we have learned from our calculations of 
screening in a hot wind. We compare our calculation of the jet quenching parameter q to the very 
different approach to energy loss in Refs. ^3), ^ |25], ^6], |2?[ in Section 6.2. In Section 6.3, we 



compare our calculation of q in M = 4 SYM to that extracted from RHIC data. Given that we find 
surprisingly good agreement between qsYM and that extracted from RHIC data, in Section 6.4 we 
enumerate the differences and similarities between M = 4 SYM and QCD. Finally, in Section 6.5 we 
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collect what is known about how q changes from the quark-gluon plasma of one gauge theory to that 
of another, including deriving a new result which allows the determination of q in any conformal 
theory with a gravity dual. We use this result to estimate (}qcd I qsYM- 

6.1 Velocity dependence of screening length 

We can summarize what we have learned from our calculation in Sections |3] and |5] of the quark- 
antiquark potential and screening length in a hot wind as follows. We find that the screening length 
L miD[ of an external quark in an M = 4 SYM plasma with velocity v = tanh rj can be written as 

£max = TfTT \f= > (6-1) 

ixl vcoshr? 

where 9 is the angle between the orientation of the dipole and the velocity of the moving thermal 
medium in the rest frame of the dipole. f(rj, 9) is only weakly dependent on both of its arguments. 
That is, it is close to constant. In fact, for any values of rj and 9, f(rj,6) lies between 0.74 and 
0.87. The limiting cases are f(rj = 0) ~ 0.87 for all 9, and /(oo, |) ~ 0.74. For a given rj, f(rj,9) 
is a monotonically decreasing function of 9 as 9 varies from to tt/2. For a given 9, f(r),9) is a 
monotonically decreasing function of rj. As rj — > oo, we find f(r},9) = /(oo,#)(l + 0(1/ cosh 2 r/)). 

For M = 4 SYM theory, since the energy density e oc T 4 , in the large rj limit equation (|6.1| ) can 
also be thought of as 

L ma , x oc i , (6-2) 
£(77)4 

where e(rj) = cosh 2 7?£:(0) is the energy density of the boosted medium. 

As discussed in Ref. H31H , if the velocity scaling of L max that we have found, namely ( |6.1D and 
( B.2|) , holds for QCD, it will have qualitative consequences for quarkonium suppression in heavy 
ion collisions at RHIC and LHC Since our discussion in earlier sections only involves the AdSs 
part of the geometry, the scaling fl6.2| ) applies to any conformal field theory with a gravity dual at 
finite temperature. To the extent that the QGP of QCD at RHIC temperature is close to being 
conformal, one is tempted to view this as a support of the applicability of (|6.1|) and ( |6.2|) to QCD. 
The results of Ref. [|30] further support this view. These authors studied large-spin mesons in a hot 
wind in a confining, nonsupersymmetric theory and found that they dissociate beyond a maximum 
wind velocity. The relation between the size L of these mesons and their dissociation velocity v is 
consistent with L oc (1 — v 2 ) 1 ^. 

For more general theories with a gravity dual, one can use the generic metric ( |6.27| ) which we 
introduce below to study the screening length. A nice argument presented by Caceres, Natsuume 
and Okamura in Ref. []33] indicates that in the large 77 limit, one would generically have 

imax OC — ^— (6.3) 

for some index v. In particular, for any gauge theory which is dual to an asymptotically AdSs 
geometry, one would find v = j as in ( |6.2| ). Examples include M = 4 SYM with nonzero R- 



charge chemical potentials, studied in Refs. [33, (Note that since chemical potentials introduce 
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additional mass scales, the dependence of L mSLX on temperature is rather complicated and it is (|6.2|) 
which generalizes, not (|6.1p.) 

For non-conformal theories the scaling index v can deviate from 1/4. One measure of the 
deviation from conformality is the deviation of the sound velocity from the conformal value of 
1 / y/3. Following a similar argument in Ref. [lfj concerning the value of q in non-conformal theories, 
Caceres, Natsuume and Okamura suggested that for theories which are close to being conformal, 
the index v may profitably be written as 



1 

v = - + c5 + 



1 

3 



(6.4) 



with c some constant. For the cascading gauge theories of Ref. |p§ 

,2 



1.5T C [§S 



c = —9/16 meaning that if 
the index is v ~ 0.21-0.24. 



Vg ~ 0.27—0.31 in these theories as is the case in QCD at T 
This suggests that for QCD the scaling (|6.2j ) with v = 1/4 should be a very good guide. 

As discussed at the end of Section 3.1, the velocity scaling ( p.2\j describes how screening lengths 
and, correspondingly, dissociation temperatures drop for quarkonia moving through the thermal 
medium with some relative velocity, and so should be included in the modelling of quarkonium 
suppression. The pr-dependent pattern of quarkonium suppression predicted by ( |6.2| ) will be tested 
in future heavy ion experiments at RHIC and the LHC. 



6.2 Comparison with energy loss via drag 

Before comparing our result for q with that extracted from RHIC data, which we shall do in Section 
6.3, we discuss the differences between our approach and recent calculations |23|, |24|, |2B], ^7 
of how an external quark loses energy while being dragged through an Af = 4 SYM plasma. In 
Appendix A (see (|A.9| )) we have reproduced the "dragging string solution" of Refs. [p3 , 24 1 for 
a string attached to and trailing behind an isolated moving test quark, since we need it as the 
subtraction term in our calculation of the quark- ant iquark potential in Section 3.1. 

To describe the results of Refs. p3l |24|, p| in their proper context, let us start with the rela- 



tivistic generalization of the Langevin equations for a quark moving through some thermal medium 
(see for example Ref. p9|, |2"5f) 



dp L 
dt 

dp T 
dt 



frit) , 



(6.5) 
(6.6) 



where pi and pt are the longitudinal and transverse momentum of the quark, respectively. (We 
have simplified the notation by dropping the spatial indices on transverse quantities.) Henceforth 
we shall denote pl by p. fa and fa are random fluctuating forces in the longitudinal and transverse 
directions, which satisfy 



(fa(t)fa(t')) = K L (p)6(t-t') 

^T(t)fa(t'))=K T ( P )6(t-t') 



(6.7) 
(6.8) 
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Kiip) and two times kt{p) describe how much longitudinal and transverse momentum squared is 
transferred to the quark per unit time. Note that at zero velocity, Kl(0) = «t(0) whereas for p > 
one expects that kl(p) 7^ k t(p)- Also, upon assuming that the momentum fluctuations of the 
particle are in equilibrium with the thermal medium, as appropriate at zero velocity, a fluctuation- 
dissipation theorem relates p,(0) to /«l(0) via the Einstein relation 

where m is the static mass of the quark. The relation (|6.9|) is also not expected to hold for p > 0. 



(See Ref. |99| for examples. 



The central result of Refs. [E3|, is that an isolated quark moving through the M = 4 SYM 



plasma with its trailing, dragging, string feels a drag force proportional to p, described by a mo- 
mentum independent drag coefficient 

Independently, kl(0) in the M = 4 SYM plasma was calculated directly in Ref. [2l|. These authors 



found that kl(0) and (i(0) in (|6.10| ) indeed satisfy ( |6.9|) , which can be considered a new consistency 



check of the AdS / CFT framework. 

Eq. (|6.10D demonstrates that in the high energy limit the energy loss mechanism in strongly 
coupled M = 4 SYM theory is very different from that in QCD. The origin of the difference is the fact 
that Af = 4 SYM is not asymptotically free. In QCD, as discussed in Section ||, the average energy 
loss of a parton in the high energy limit is independent of p (or at most logarithmically dependent 
on p) and it is proportional to the square of the distance travelled through the medium. There, the 
dominant mechanism by which a high energy parton loses its energy is through radiating gluons 
which have a high enough transverse momentum kr (and an even higher energy) that a s evaluated 
at kx is weak, meaning that the dominant energy loss processes can be described perturbatively, 
with nonperturbative physics at scales of order the temperature coming in only via the description 
of the repeated soft interactions between the radiated gluon and the medium, and between the 
original high energy parton and the medium. The effect of these nonperturbative soft interactions 
is encoded in the jet quenching parameter q, which can be defined nonperturbatively via a light-like 
Wilson loop as described in Section 2. 



In Ref. |15j and in the present paper, we seek insights about q in QCD by calculating the 
analogous quantity in Af = 4 SYM. We do not attempt to describe the full process of energy loss in 
Af = 4 SYM because the asymptotic freedom of QCD is crucial to the description of the radiative 
parton energy loss process which dominates at high energies, making it impossible to model the 
physics of QCD parton energy loss at high energies in a theory like Af = 4 SYM which is strongly 
interacting at all scales. 

Even though the drag coefficient fl6.10| ) describes energy loss in a A" = 4 SYM plasma even for 



quarks moving relativistically ||23|| , it cannot be used to extract q or k,t{p oo). As we remarked 



after equation (|6.8|), except in the low- velocity limit one does not expect nrip) to be equal to kl(p) 
or kl(p) to be related to fi(p) via the Einstein relation. Indeed, the direct calculation of k,l in 
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Ref. [25] manifestly requires modification at nonzero p [ |f 00|| . The «l(0) found in Ref. |25| and 



indirectly in Ref. |23[] has the same parametric dependence on A and T as q in ( |6.11D , and is smaller 
than q by a purely numerical factor of ~ 1.20. This is curious, since the quantity kl(0) = kt(0) has 
no evident relation to jet quenching (since jets by definition are relativistic) or to the jet quenching 
parameter q. 

To look for connections between M = 4 SYM energy loss described by ( |6.1(J| ) and data from 
RHIC one should then seek circumstances in which all aspects of the energy loss process are strongly 
coupled. Perhaps the energy loss of quarks which are slowly moving and yet energetic — i.e. quarks 



which are heavy — is the best example, as stressed by many of the recent papers |23], |24|, |25|, g7) , 
although extracting the contribution from energy loss to the medium modification of charmed 
meson production in the regime in which the parent c-quark is slowly moving presents considerable 
challenges. Furthermore, precisely because the entire energy loss process is treated in this approach, 
further questions like how the M = 4 SYM medium responds to the dragging quark (i.e. "where 



does the lost energy go?") can be addressed |27 



6.3 Comparison of q of M = 4 SYM with experimental estimate 

We turn now to the comparison between q in M = 4 SYM with that extracted from current RHIC 
data. In Eq. ( |3.45| ) we found that 



7r 3/2 r (-) 

Qsym = ( ^ 4J v^T 3 « 26.69 VasruNcT 3 . (6.11) 

Taking N = 3 and «sym = o> thinking ojqcd = \ reasonable for temperatures not far above the 
QCD phase transition, we shall use A = 6n to make estimates. 5 From ( |6.11|) , we find 



Qsym = 4.5, 10.6, 20.7 GeV 2 /fm for T = 300, 400, 500 MeV . (6.12) 



This qualitative comparison will suffice for our purposes. The question of how to relate the coupling strengths 
of different thermal quantum field theories to each other is not guaranteed to have an unambiguous answer. Here 
we note that, unlike asYM, the value of cxqcb(T) is not well-defined in a strongly interacting quark-gluon plasma. 
It may thus be preferable to make this comparison using a quantity calculated nonperturbatively in both theories, 
and compare the QCD value for the QGP at a few times the QCD critical temperature T c , without reference to any 
ttQCD, to the A-dependent SYM value, thus fixing A. Because we want to keep the question of how to fix A separate 
from the question of how the difference in the number of degrees of freedom in QCD and TV = 4 SYM affects q, it is 
also important to choose a quantity which is independent of the number of degrees of freedom. (For our purposes 
this rules out using the Debye mass, as suggested in Ref. |x| in a different context.) One recent proposal is 
to compare the shapes of the static quark-antiquark potential. The screening length L max is independent of A in 
JV = 4 SYM, so that cannot be used. Instead, one has to compare the shapes of the potentials themselves, which is 
not straightforward. Another possibility is to use the lattice QCD calculation of the ratio of the energy density of 
the QCD quark-gluon plasma to that of a non-interacting Stefan-Boltzmann gas of quarks and gluons. According to 
lattice QCD calculations done with two or three flavors of quarks, this ratio rises to about 0.78-0.82 for T ~ (1.5 — 2)T C 
and then flattens at higher T |J2). In SYM, it is | in the A — > oo limit |)3) and, using the leading correction to 
this result which is +||C(3)A -3 / 2 the range 0.78-0.82 corresponds to 9 < A < 15. The pressure in two- and 
three-flavor QCD also approaches about 0.8 times its Stefan-Boltzmann value, but only at the higher temperature 
T ~ (2.5 — 3)T C [p2|, suggesting that a comparison of this sort could be more appropriate at the LHC than at RHIC. 
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We have computed ( 6.11 ) in the large N and large 't Hooft coupling A limit and thus it is 
only the leading order term in a double expansion in 1/y/X and 1/N 2 . There are two sources of 
contributions to 1/vA corrections: from the fluctuation of the string world sheet (|3.5| ) and from 



the modification of the background geometry (3.8) due to a' corrections. The authors of Ref. 



find that the leading contribution of the second type of correction is given by 

qsYM{X) = qsYM (l - 1.765A~i + • ■ •) . (6.13) 

If one uses «sym = \, then the correction from ( |6.13| ) is about 2%. The corrections due to fluc- 
tuations of the world sheet, the leading order term of which is expected to be of order l/\f\ ) are 
harder to compute and are not known at the moment. 1/N 2 corrections, which require string loop 
calculations, are also beyond currently available technology. It is worth noting, however, that the 
ratio of energy density and pressure to their Stefan-Boltzmann values are both quite insensitive to 
changing N from 3 to 4 to 8 [95]. 

In a heavy ion collision, q decreases with time r as the hot fluid expands and cools. The time- 
averaged q which has been determined in comparison with RHIC data is q = jj^2 f^ +L ^ T <?( r ) dr, 
found to be around 5-15 GeV 2 /fm [^, 97|. If we assume a one-dimensional Bjorken expansion with 

1 /3 

T(t) = T (^) , take tq = 0.5 fm, and take L~ / v^2 = 2 fm, the estimated mean distance travelled 



in the medium by those hard partons which "escape" and are detected H97|| , we find that to obtain 



q = 5 GeV 2 /fm from (|6.11|) we need To such that T(l fm) « 310 MeV, only slightly higher than 
that expected from hydrodynamic modelling |98[] . With L~ /y2 = 1.5 fm, we find q = 5 GeV 2 /fm 
for T(l fm) m 280 MeV, fully consistent with expectations. There are currently too many uncer- 
tainties in the various components of this comparison to make a strong statement, but it seems 
clear that the q given by ( p.ll| ) that we have calculated in the quark-gluon plasma of M = 4 SYM 
is in qualitative agreement with q = 5 GeV 2 /fm, which is in turn consistent with RHIC data. 

The value obtained from ( |6.11| ) assumes that the medium is static. However, in a relativistic 
heavy ion collision, the medium itself develops strong collective flow, meaning that the hard parton 
is traversing a moving medium — it feels a wind. Thus to compare (|6.11| ) with the experimental 



estimate we should include the effects of the wind on q that we discussed in Section [|. We found 
in Eq. [113] that 

q = 7/(1 - v f cos9) q , (6.14) 

where Vf = tanh?7j is the velocity of the wind, 7/ = 1/ ^1 — v 2 , and 6 is the angle between the 

direction of motion of the hard parton and the direction of the wind, go is the value of q in the 
absence of a wind. The result ( |6.14| ) for the dependence of q on collective flow is valid in QCD 
and in M = 4 SYM and in the quark-gluon plasma of any other gauge theory, since its derivation 
(see Section 4) relies only on properties of Lorentz transformations. If we crudely guess that head 
winds are as likely as tail winds, and that the typical transverse wind velocity seen by a high energy 
parton is about half the speed of light, q is increased relative to that in ( |6.11| ) by a factor of 1.16. A 
credible evaluation of the consequences of (|6.14|) for the time-averaged q extracted from data will, 



however, require careful modelling of the geometry of the collision and the time-development of the 



collective flow velocity, as in Refs. [85. 86. 
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The weak-coupling QCD estimate of q given in ( 2.18 ) is, when evaluated with a s = 1/2, smaller 
than that in fl6.12| ) by about a factor of 5. This means that in order to make the weak-coupling 
estimate consistent with data from RHIC, we would have to choose a s > 1, certainly beyond weak- 
coupling. The strong coupling calculation of q in Af = 4 SYM is certainly in better agreement with 
data from RHIC than the weak-coupling calculation in QCD. The obvious question, then, is how 
much the strong coupling result will change as one modifies the theory from Af = 4 SYM towards, 
and ultimately to, QCD. We shall describe the current state of our ability to answer this question 
in Section 6.5. 

The BDMPS description of radiative parton energy loss, with the nonperturbative physics of the 
medium entering through the jet quenching parameter q, is appropriate in the high parton energy 
limit. It is only by comparison to data that we can learn whether the jets being quenched at RHIC 
are sufficiently energetic for their energy loss to be described well by this formalism. To date this 
comparison has been broadly successful, albeit with q seen as a free parameter. If we understood 
the QCD prediction for q as a function of temperature in a strongly interacting quark-gluon plasma 
even at the factor of two level, this would make this comparison more stringent. And, it would 
turn q into a "thermometer" with a calibration error of only a factor of 2 1 / 3 , which would be an 
exceptionally valuable addition since one of the biggest current weaknesses in our understanding of 
RHIC phenomena is that we do not have an experimental measure of the temperature at, say, a time 
of 1 fm after the collision. In the next two subsections, we first frame the questions that need to be 
thought through if it is to be possible to go from our calculation of q in Af = 4 SYM to a factor of 
two understanding of q in QCD, and then review and extend the calculations of q in various other 
gauge theories, yielding a conjecture for how to estimate q in QCD at a semi-quantitative level, still 
with many caveats. Our conjecture is that (Jqc d / qsY m is of the order of -y/47.5/120 ~ 0.63, namely 
the square root of the ratio of the numbers of degrees of freedom in the two theories. 

6.4 Af = 4 SYM versus QCD 

We found in Section 6.3 that q calculated in Af = 4 SYM theory is close to the value extracted 
from RHIC data. Given that RHIC is probing the quark-gluon plasma of QCD, is this agreement 
meaningful or accidental? In what respects can the strongly interacting plasma of Af = 4 SYM 
theory give a reasonable description of the quark-gluon plasma in QCD? After all, at a microscopic 
level Af = 4 SYM is very different from QCD: 

• The theory is conformal, supersymmetric and contains additional global symmetry. The 
coupling does not run and there is no confinement. 

• No dynamic quarks, no chiral symmetry and no chiral symmetry breaking. 

• Additional scalar and fermionic fields in the adjoint representation. 

These features of course make the vacuum sectors of the two theories very different. However, if 
the quark-gluon plasma in QCD is strongly interacting, as indicated by data from RHIC, then one 
may ask whether the macroscopic properties, both thermodynamic and dynamic, of quark-gluon 
plasma at sufficiently strong coupling may be insensitive to differences between the theories which 
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seem stark in vacuum. It is often the case that macroscopic properties of a sufficiently excited 
many-body system are not sensitive to the detailed underlying dynamics, with systems within the 
same universality class exhibiting similar phenomena. We are used to the idea that all metals, or all 
liquids, or all ferromagnets have common, defining, characteristics even though they may differ very 
significantly at a microscopic level. What we are asking here is what are the defining commonalities 
of quark-gluon plasmas in different theories, and in what instances do these commonalities allow 
qualitative or semi- quantitative lessons learned about the quark-gluon plasma of one theory to be 
applied to that of another. 

Returning to the differences between QCD and M = 4 SYM, many are obviously irrelevant to 
a comparison between strongly-coupled plasmas in the two theories. After all, supersymmetry is 
explicitly and badly broken at high temperature and, above T c in QCD, there is no confinement 
and no chiral condensate. Furthermore, in a strongly interacting liquid there are, by definition, 
no well-defined, long-lived quasiparticles anyway, making it plausible that observables or ratios 
of observables can be found which are insensitive to the differences between microscopic degrees 
of freedom and interactions. Because Af = 4 SYM is a conformal theory whereas QCD is not, 
Af = 4 SYM cannot be used to describe QCD at or below its phase transition at T ~ T c , and 
cannot be used to describe QCD at asymptotically high temperatures. However, there are a variety 
of indications from lattice QCD calculations (enumerated below) that QCD thermodynamics is 
reasonably well approximated as conformal in a range of temperatures from about 2T C up to some 
higher temperature not currently determined. It is not currently known whether the quark-gluon 
plasma of QCD, as explored at RHIC and in lattice QCD calculations, and that of Af = 4 SYM, as 
explored using AdS / CFT calculations, are in the same universality class, or even in what sense this 
question could be made precise. However, given the rapidly increasing list of similarities between 
the two quark-gluon plasmas, it does not seem too far-fetched to imagine. Here is a list of some 
of the similarities between the quark-gluon plasmas of the two theories, notwithstanding the stark 
differences between their vacua, with thermodynamic comparisons listed first followed by dynamic 
comparisons: 



Above about 1.2T C , the ratio of the energy density e in 2- and 3-flavor QCD to that in the 
absence of interactions is close to T-independent and takes on the value of about 0.8 [|101|, j92|, 



89 

the A 



In zero-flavor QCD, this ratio is closer to 0.9 [ |102|| . In Af 
oo limit |3[ and is 0.8 for A ~ 11 @]. 



4 SYM, this ratio is 3/4 in 



Above about 2.5T C , the ratio of the pressure P in 2- or 3-flavor QCD to that in the absence of 



In Af 



interactions is also about 0.8. In zero-flavor QCD, this ratio is closer to 0.9 ||102 
SYM, this ratio must be the same as that defined via the energy density, a condition that 
is satisfied well in QCD. Note, however, that for T < 2.5T C , the deviation from conformality 
parametrized by e — 3P is significant. This suggests that the use of a conformal theory like 
Af = 4 SYM as a model for the quark-gluon plasma may be more quantitatively reliable for 
heavy ion collisions at the LHC than at RHIC, since RHIC is likely exploring temperatures 
that are less than 2T C whereas the LHC can be expected to reach temperatures that are higher 
by about a factor of two. 
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All the results in M = 4 SYM have been obtained in the N —>■ oo limit. Although corrections 
are expected to be of order l/N 2 , they have not been computed. It is therefore very useful to 
test how the quantities on the QCD side of these comparisons change with N. One example 
of such a test is the calculation of Ref. ||95|| , which finds that the ratio of the pressure to its 



noninteracting value changes very little as N is changed from 3 to 4 to 8. 

The square of the speed of sound in the QCD quark gluon plasma is close to 1/3, the value for a 
conformal theory for T > 2T C ; at T = 1.5T C , it is already ~ 0.27-0.31 [[102], [103], [104], [105], |9 



The screening length defined by the potential between a quark and antiquark at rest is 
0.869/ttT ~ 0.28/T in M = 4 SYM in the large N and A limits, as calculated in Section 
3 and as calculated first in Refs. P5[ . In QCD the screening length is not sharply defined, 
since the potential does not change suddenly to zero, but operational definitions exist in the 



literature. For QCD with zero | 106 | and two [|107]| flavors, it is ~ 0.7/T and ~ 0.5/T, re- 
spectively. QCD and M = 4 SYM are therefore qualitatively comparable in this regard, with 
the quantitative difference between them plausibly reflecting the larger number of degrees of 
freedom in Af = 4 SYM. 



Turning now to dynamic quantities, the shear viscosity in units of the entropy density is 1/47T 
in J\f = 4 SYM in the A — > oo limit ||, and is ~ 1.25/47T for A = 6ir @. Given the degree to 
which data on the azimuthal anisotropy of RHIC collisions are well-described by zero- viscosity 
hydrodynamics, the ratio of the shear viscosity to the entropy density has been estimated to 



be comparably small in the quark-gluon plasma at RHIC [ll|. A quantitative extraction 
of 7] from RHIC data requires viscous hydrodynamic calculations, which are currently being 
pursued by various groups | |108| . 

In M = 4 SYM, the quark- antiquark screening length scales with velocity according to L s (v) ~ 
-^s(O) / a/7, as discussed in Ref. |3T[| and Sections 3.1 and 5 above. Comparison of this predicted 



scaling to QCD awaits data from RHIC on the py-dependence of J/\l/ suppression at RHIC 
at pt > 5 GeV and on the pattern of py-dependence of J/\l/ and Upsilon suppression at the 
LHC. 

As we have discussed in Section 6.3, the jet quenching parameter q in M = 4 SYM is close to 



the value extracted from RHIC data 15 



To understand whether the above similarities are meaningful, one avenue is to study strongly 
interacting quark-gluon plasmas in other non-Abelian gauge theories with dual gravity descriptions 
and see whether a general picture emerges. In the case of the ratio of the shear viscosity to the 
entropy density, it was indeed found that this ratio is the same in a broad class of gauge theories j|, [|. 
A necessary (but not sufficient) condition for this striking lack of dependence on the nature of the 
microscopic theory is that the dependence on the number of degrees of freedom cancels in the 
dimensionless ratio. In thinking about how the value of q, a dimensionful quantity, may change 
in going from M = 4 SYM to QCD, it seems to us that the two most pressing questions are how 
q depends on the number of degrees of freedom and on the fact that QCD includes fundamentals 
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whereas all the degrees of freedom in M = 4 SYM are adjoints. These seem to have greater 
potential to change q significantly than do 1/y/X corrections, 1/iV 2 corrections or corrections due 
to the deviation away from conformality. We cannot currently address the effect of fundamentals. 
In Section 6.5, we consider how q depends on the number of degrees of freedom. 

6.5 The jet quenching parameter and degrees of freedom 

One qualitative feature of Eq. ( |6.11|) is that at strong coupling q is proportional to y/X, not to the 
number of degrees of freedom ~ iV 2 . This means that at strong coupling, q cannot be thought of 
as "measuring" either the entropy density s or what is sometimes described as a "gluon number 
density" or e 3 / 4 as had been expected [[F(| , since both s and the energy density e are proportional to 
N 2 X°. Whereas the ratio of the shear viscosity to s turns out to be universal in theories with gravity 
duals, the ratio q/s vanishes in the large- N limit. Even though (|6.11|) upends prior intuition on this 
point, it nevertheless seems that q should have some straightforward dependence on the number of 
degrees of freedom in the theory. We shall now try to make this intuition precise. 

We first examine how q of different conformal field theories with a type IIB supergravity dual 
compare to each other, using that of Af = 4 SYM as a reference point. The ten-dimensonal metric 
dual to a conformal field theory at zero temperature can be written in the form 

dsj = n 2 (y)R 2 (ds 2 AdS5 + ds 2 Ms (y)) , (6.15) 

where R is the curvature radius of AdSs and the metric for AdSs inside the parenthesis is normalized 
to have curvature radius unity. ds 2 Ms is the metric of an internal five-dimensional manifold. The 
warp factor Q 2 (y) depends only on the coordinates y of the internal manifold. 6 To put the theory at 
finite temperature, one replaces ds 2 AdSb by the metric of ( |3.8| ) of an AdS-Schwarzschild black hole. 
The computation of q is identical to what we have done before and we find that (cf. ( |3.45|) ) 



7T 



3 / 2 r (f) Q 2 (y)R 



■2 



r(!) 



a' 



T . (6.16) 



We shall compare different theories at fixed values of N, X = A7ig s N and T. Note that the M = 4 
SYM relation vA = R 2 jol is modified in this more general context, as we shall see below in (|6.17f) . 



Let us first consider theories in which the warp factor is trivial, i.e. Q = 1. In addition to £5, 
which corresponds to M = 4 SYM theory, an infinite number of examples of such dual pairs are 
now known |109| , |110(| , with the boundary conformal field theories being quiver gauge theories with 



product gauge groups. In the simplest example, the manifold M5 is a manifold with the topology 
of £5 known as T 1 ' 1 and the corresponding boundary theory is the Klebanov-Witten CFT [|111| . 



In general, Type IIB supergravity equations of motion fix the curvature radius R in terms of the 
number iV of D3-branes as 

R 4 = 4irg s Na' 2 ^- , (6.17) 



6 The full supergravity solution also involves a self-dual five-form, which we will normalize to have the flux of N 
D3-branes, and possibly three- forms (when f2 is nontrivial). 

7 In this case dsj^ is an Einstein manifold with curvature RqP — 45^ . If one requires the boundary theory be 
supersymmctric, then M$ needs to be a Sasaki-Einstein manifold. 
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where u>m 5 and Ug 5 = tt 3 are the volume of ds\ Is and ds 2 s respectively. Plugging (]6.17|) into ( |6.16| ) 
we find that (with Q = 1) 



Qcft 



us 5 



Noting that the central charge a of a CFT can be written as ||112| , |1 13 

acFT 



we can rewrite ( |6.18|) as 



Qcft 



us 5 

W M 5 



ClcFT 



qM=A V a A/"=4 

For example, in the Klebanov-Witten theory | |111| , 



(6.18) 



(6.19) 



(6.20) 



Qkw 
QJV=4 




(6.21) 



As another example, note that if you start with M = 4 SYM, described by Type IIB string theory 
on AdSs x 5*5, and orbifold the S$ by Z 2 , the central charge of the CFT doubles and q increases by 
a factor of \/2. To understand the implications of ( |6.18 ), recall that the entropy density of a CFT 
is related to its central charge such that 



SCFT _ &CFT 
SjV=4 a Af=4 



(6.22) 



making it clear that the central charge counts the number of thermodynamic degrees of freedom in 
the theory. We conclude that in any conformal theory with a gravity dual ( |6.15| ) with Q = 1, 



Qcft 
QJV=4 



SCFT 



(6.23) 



Note that even though q oc VX'V and s oc iV 2 A , these factors cancel in the ratios on the left and 
right hand sides of (|6.23|). This equation should be read as saying that, in the relevant class of 



theories, q/VX is proportional to \Js/N 2 . 

When the warp factor Q(y) in (|6.15|) is nontrivial, the value of q depends on where in the 
internal space M5 we put the probe brane. In these theories, different types of quarks, which 
correspond to putting branes at different locations in the internal manifold, have different values of 
q, as was first pointed out in Ref. [fT7|| . An example of ( |6.15| ) with nontrivial Q is the Pilch- Warner 
geometry ||114|| which is dual to the = 1 superconformal CFT of Leigh and Strassler ||115|| . In 
this case the computation of the normalization condition and central charge are more complicated 



(see, e.g., Ref. 



but the result is again simple to state: 



qLs{y) 
qN=A 



(6.24) 
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with the ratio of <fs again proportional to the square root of the ratio of central charges. For the 
Pilch- Warner geometry, 8 = §| and 

s /3 — COS 26 \ 5 7T . . 

fi 2 (y)=(^ g J' °- e -2' (6 ' 25) 

where # is an angle specifying a location within the internal M 5 . 

Thus for any two CFTs with a holographic dual, the ratio of their jet quenching parameters q 
is proportional to the square root of the ratio of their central charges, and hence to the square root 
of the ratio of their number of degrees of freedom. "Proportional to" becomes "equal to" if Q = 1. 
In particular, if two CFTs are connected by a renormalization group flow then q for the UV theory 
is always larger than that of the IR theory. 

Since QCD is not a CFT, we cannot directly apply the result we have just derived to QCD. 
However, to the extent that the quark-gluon plasma of QCD is approximately conformal, as we 
have discussed in Section 6.4, perhaps our result for comparing CFTs can be used as a guide. In 
doing so it seems fair to set — 1 since there is no indication that if QCD with Nf massless flavors 
of quarks had a holographic dual, there would be quarks with differing physics corresponding to 
branes at differing locations in an internal manifold. So, we conjecture that 



Qqcd _ /^_./ 47 - 5 ~o.63 (6.26) 



QM=4 V s a^=4 V 120 

is a good estimate of the effect of the difference between the number of degrees of freedom in the 
two theories on q. We have used N = 3 in both theories, and have used Nf = 3 in QCD. 9 And, 
we make the comparison with A chosen in the Af = 4 SYM theory such that the ratio of s to its 
value in a noninteracting theory is the same as that in the QCD quark-gluon plasma. It would 
be good to ask how ( |6.26| ) is affected by the fact that some of the degrees of freedom in QCD are 
fundamentals. Unfortunately, we do not currently have any examples of calculations of q in theories 
with fundamentals among the degrees of freedom of the strongly interacting plasma. 

Next, we ask how q is affected by deviations from conformality. In a nonconformal theory with 
a dual gravity description, the bulk metric in the string frame can generically be written in the form 

dv^ 

ds 2 = g(r, y) [-(1 - f(r, y))dt 2 + dx 2 )] + —— - + ds\ h (y, r) , (6.27) 
where y again denotes coordinates of the internal manifold. The corresponding q can be written 

sin 



as 




9 = —, / • ( 6 - 28 ) 



Note that our normalization ( 6.1 5| ) is different from that in | 114 | and 116 ] 



9 In QCD, the gluons contribute 2(N 2 - 1) and the Nf flavors of quarks contribute ^ANNf. In Af = 4 SYM, the 
gauge bosons, fermions, and scalars contribute 2{N 2 — 1), |8(iV 2 — 1) and 6(iV 2 — 1), respectively. 
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It appears hard to extract a general story from Q6.28 ) directly, without studying more examples. 
We can write q in the form 

q=^a(X,^)^T 3 (6.29) 

with jj, some mass parameter(s) of the nonconformal theory. Motivated from our discussion of 
conformal theories it is tempting to speculate that a(A, ^) can be considered as a measure of the 
number of degrees of freedom of a theory at an energy scale T. The jet quenching parameter q 
for the nonconformal cascading gauge theories of Ref. []88| is known in the high temperature limit 
and the result is consistent with the hypothesis that a(A, ^) decreases with renormalization group 
flow [jl6 1 . It would also be interesting to compute q for the geometry discussed in [|1 1 7| to see 



whether the function a(A, ^) decreases monotonically with renormalization group flow. In other 
words, the function a(A, ^) defined by the jet quenching parameter is a candidate resolution of the 
long-standing challenge to find a four dimensional analogue of the c— function of two dimensional 
conformal field theory. 

We can also ask seek to evaluate how much q is affected if the theory is "as nonconformal" as 
the quark-gluon plasma of QCD is at a few times its T c . There is no one prescription for quantifying 
nonconformality. However, the analysis of the cascading gauge theories of Ref. |8~8l provides a nice 



example, as the effect of the nonconformality on q can be written |T6" | 



^cascading o i n / 2 



Qkw V \3 



3.12 --^ , (6.30) 



where v s is the speed of sound. For v% in the range 0.27-0.31, as in QCD at T = 1.5T C |jS9 , the 
effect of the deviation from nonconformality on q ranges from 6% to 18%. 

We have obtained one nontrivial check in a nonconformal theory of our conjecture that a, 
defined from the jet quenching parameter via Q6.29| ), is a measure of the number of degrees of 



freedom, as required if our specific conjecture ( |6.26| ) is to hold. Consider q for (p + l)-dimensional 



super- Yang-Mills theories (with 16 supercharges) living at the boundary of the geometry describing 
a large number of non-extremal black Dp-branes ||119|| . We will restrict to p < 5. The case p = 3 



is M = 4 SYM; the cases p = 2 and p = 4 correspond to nonconformal theories in 2 + 1- and 
4 + 1-dimensions. The metric dual to these theories can be written as 

ds 2 = JA^pIl (^fdf + ds l + (_?J) 2 ^ + z 2 dnl p ) , (6.31) 

where 10 



z 2 \ \ 5 — p J f 



A = AN, /=l-(f) • d r = 2'- 2 ^r (IZ£) . (6.32) 



10 Here we are following standard string theory convention and normalizing the gauge coupling constant as g\j = 

3-p 

{2-k) p g s a' 2 . For p = 3, the gauge coupling in the standard field theoretical convention which we have used 



9ym = ^9m' meaning that A = |A. For p = 3, the relation between these coordinates and those in (3.8) 
is z = R 2 /r with R 2 — a'y/X, meaning that z — R 2 /r and / = R 2 f/r 2 . 
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dsi is the metric for flat p- dimensional Euclidean space and dVL 2 d is the metric for a (^-dimensional 
sphere. Note that g 2 Ml and hence A, have mass dimension 3 — p. The horizon is at z = zq and the 
boundary is at z = 0. The temperature can be obtained as 

7 - v 1 

T = y - . (6.33) 

5 — p 2-kzq 

The energy density e and entropy density s of the systems can be written as 

e = N*±^b p X*jT (T)l^ 1 , (6.34) 



(T)T P 



(6.35) 



where 



A eff (T) = AT P ~ 3 , 



(7-p) 



7-p 



(6.36) 



Aeff(T) is the effective dimensionless coupling at temperature T. Note that equation ( |b.35| ) indicates 

P-3 

that the quantity b p \l^ v (T) characterizes the number of degrees of freedom at temperature T. By 
following the procedure of Ref. |15| or Section |3] above, or by simply applying ( |6.28|) to ( |6.31|) , we 
find that 



Q 



1 p-3 



' 5-p \ 
14- 2p^ 



6|A^-(T)V%^T)T 3 . 



(6.37) 



We see from ( 6.3?] ) that the quantity a(A e ff) defined as in OPD has the same dependence on A e fr 
that the entropy density (|6.35|) has. 

The calculation of q in the nonconformal p = 2 and p = 4 Dp-brane theories supports our con- 
jecture that a, defined from q via ( |6.29| ), measures the number of degrees of freedom at temperature 
T. 11 This conjecture can be further tested by computing q in other nonconformal theories, like for 
example the M = 2* theory of Refs. | |1 |5|. This would also allow us to test our conjecture that 



11 A relation analogous to that between q and the number of degrees of freedom may also be valid for the longitudinal 
drag coefficient (J-l oi Refs. |23|, [ 
( |6.29 ) we can introduce 



defined in (6.5) and given by (6.10) in Af = 4 SYM. In particular, following 



Mi 



m\ 1 



(6.38) 



and o(A, ^) and a(X, ^) could coincide up to some numerical constant. This holds for all the CFTs with a IIB 
supergravity dual, meaning that in all these theories a is proportional to the number of degrees of freedom just as a 
is. This is also the case for the Dp-brane, where we find that 



ML 



1 (h-p 



I - V 
5 - p 



1 P-3 

2 5-p . 



(6.39) 
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the function a(A, ^) decreases under renormalization group flow, and it would test the conclusion 
indicated by fl6.30| ) that deviations from conformality with a magnitude comparable to those in the 
QCD quark-gluon plasma do not change q much. 

If our conjecture relating q to the number of degrees of freedom, valid for any conformal theory, 
survives being further tested via the calculation of q in more examples of nonconformal theories 
with gravity duals and, even better, in theories with fundamentals we will then have a new example 
of a common feature of strongly interacting quark-gluon plasmas. Furthermore, the conjecture 
( |6.26| ) together with our result ( |6.11| ) for (}sym will then provide a theoretical prediction for the jet 
quenching parameter q in the strongly interacting quark-gluon plasma of QCD. In order to make 
quantitative contact with data, we will further need to model the effects of collective flow on q, 
using the result (|6.14 ) which is valid in QCD. On the experimental front, we can look forward 
to studies of jet quenching being extended to higher and higher transverse momentum jets as 
RHIC runs at higher luminosities and as the LHC comes on line. Furthermore, particularly at the 
LHC new observables sensitive to parton energy loss will be developed. Making these assumptions 
about (near) future theoretical and experimental developments, we can look forward to a stringent 
comparison between experimental and theoretical determinations of the jet quenching parameter in 
the strongly interacting quark-gluon plasma of QCD. 



Note added 



Shortly after the completion of our work, two papers ||120| , |121|| appeared which calculated the mean 
squared momentum transfer kt(p) in ( |6.8|) for a heavy external quark of mass M moving with a 
velocity v = tanh^ through the M = 4 plasma. These calculations are valid for the kinematic 
poj |T2TH 

\j cosh 7] < —f=— = A , 



regime 



\f\T 
00 only if the M 



(6.40) 



00 limit has been taken first. In 



meaning that they can be extended to rj 
the regime ( |6.40| ), one finds 

kt — \Zcosh~7y \/~\ it T 3 . (6.41) 

In contrast, as discussed in Section |3.3| our calculation of q requires \/ cosh rj > A. Since this 
kinematic region does not overlap with ( |6.40 ), a direct comparison of q and k,t is difficult. However, 
we agree with Refs. |120| , |121|| that it would be desirable to have a better understanding of whether 
there is a connection between q and Ky, and what this connection could be. To illustrate this open 
issue, let us make the following remarks. On the one hand, it is obvious that k,t cannot play the role 
of q within the BDMPS energy loss formalism ( |2.14|) , on which the definition of q is based. To see 
this, note first that q in ( 2.14 ) is defined for a quark moving strictly along the light-cone, rj —>■ 00, 
for which in (|6.41|) diverges. Then, note that replacing q in (|2.14j ) by a divergent quantity leads 
to an ill-defined expression. On the other hand, the value of q in ( |2.14| ) is known to set the scale 
of the transverse momentum broadening of the medium-modified gluon radiation. (For example, 
see Fig. 1 of Ref. f72|| .) So, while the value of q is known to set the transverse momentum scale of 
BDMPS energy loss, a direct connection between q and may not be straightforward. 
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We also stress that in the BDMPS formalism, the quantity q which parameterizes radiative 
energy loss of heavy quarks moving at high velocity is precisely the same as that for massless 
quarks, and is defined in terms of the short distance behavior of a strictly light-like Wilson loop 
( 12.15| ). One way of seeing this is to recall from our discussion in Section 2 that the two adjoint 
light-like Wilson lines in the Wilson loop which defines q can be thought of (loosely) as representing 
the radiated gluon (in the amplitude and in the conjugate amplitude in the calculation of gluon 
emission) which is of course on the light-cone independent of the mass of the quark. (See Section 
2 and references therein for a more complete description.) Although q itself is independent of the 
mass of the quark which is losing energy, the relation between q and AE, the average energy lost, 
does depend on the quark mass. The relation (|2.17|) , obtained from (|2.14|) , is only valid for massless 
quarks. (For the analogous expressions for massive quarks, see Ref. [ 123(1 .) The BDMPS formalism 
which relates the energy lost by massless or massive quarks to q is only valid to leading order in 
1/E. At higher order, i.e. at lower energies, the energy loss will depend on more properties of the 
medium than just the single parameter q. The divergent quantity kt computed in Refs. |120|, 1121 



cannot serve to define the jet quenching parameter q in the BDMPS radiation spectrum Q2.14 ) for 
either massless or massive quarks. 

Also after our paper appeared, Argyres et al. presented a study of "upward going" space-like 
string configurations [ |122|| . In the notation of our paper, these are found as follows. First, extend 
the AdS5 spacetime above the (no longer appropriately named) boundary D3-brane at r = Aro, 
where the quark that defines the Wilson loop is located. Then, solve (|3.32|) with the sign of y' at 
a = ±C,/2 chosen so that the solution y(a) to ( |3.32| ) begins with y(—£/2) = A, then ascends to a 
turning point at y = y m (with y m defined by (|3.33|) taking on a value which is just below \/ cosh n for 



small £ and hence small g), and then descends back down to y((/2) = A. We note that in the limit 
in which y/ cosh rj — > oo at fixed A, the turning point up to which this solution ascends is infinitely 
far above the "boundary" at A. 

these strings are probing physics at 



Hence, according to the standard IR/UV connection [79 



length scales infinitely shorter than the thickness of a Wilson line, in other words infinitely far to 
the ultraviolet of what is normally considered to be the ultraviolet cutoff in the field theory. It 
remains to be seen what field theory interpretation can be given to these upward going strings, 
but they are certainly not relevant to the evaluation of thermal expectation values of Wilson loops, 
which are located on D3-branes which bound the AdSs spacetime. 

The results of Ref. [ 122j| themselves confirm the conclusion, reached above via the IR/UV 
connection, that the upward-going strings are not related to the light-like Wilson loop which arises 
in the physics of deep inelastic scattering and radiative energy loss in QCD and which we have 
calculated in M = 4 SYM. The action of the upward-going strings of Ref. [|122|| is the same as 
that of the configuration fl3.36|) , namely 5* = za/ATLL~/2a/2, and is linear in L, the tranverse 
extent of the Wilson loop. If the expectation value of the light-like Wilson loop were to have this 
behavior, it would yield a photoabsorption probability (|2.5| ) (for the thought-experiment of deep 
inelastic scattering off quark-gluon plasma) which fell at large virtuality only like 1/Q, rather than 
the standard 1/Q 2 . [See our discussion of Eqs. ( [2.4|) and ( p75|) .1 Furthermore, if we subtract S 
from the action S, with the finite So given by ( |A.18| ) for any choice of p, this yields a negative real 
L-independent term in the exponent of (^(Cught-iikc)), making the photoabsorption probability 
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nonzero in the Q — > oo limit, meaning that a dipole of zero size would have a nonzero absorption 
probability. 
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A. Single-string drag solutions and heavy-light mesons 

The calculation of the expectation values of time-like and light-like Wilson loops in Sections 3 and 
5 required the subtraction of terms corresponding to the action of a quark and antiquark which 
propagate independently along the long sides of the Wilson loop, i.e. without seeing each other. 
In this Appendix, we enumerate the different extremal string world sheets that are possible (for 
different values of an integration constant) given a single quark at rest on the D3-brane at r = r A 
in the presence of a thermal medium moving with rapidity r\ along the ^-direction. 

If A > V cosh rj, in addition to the drag solution of Refs. |23], ^4| we find solutions in which the 



string which begins on the D3-brane at r = r A ends on a D3-brane located at r = r ^/coshr]. Such 
solutions model mesonic bound states of a heavy quark and a light quark in which the light quark 
drags behind the heavy quark. 

We discuss a single quark moving along the X3-direction. The string world sheet of this quark 
is of the form 

r = t, a = r, x 3 = x 3 (r, a) . (A.l) 

Calculating the components g a p of the induced metric (|3.6|) within this ansatz, one finds for the 
world-sheet action ( |3.5|) of two independent strings 



2 f r ° A A 2 B C 

Here, x 3 = dx 3 /dt and x' 3 = dx 3 /dr. The notational short hands / and A, B and C are defined 
in ( |3.9p and ( |3.14j ), respectively. We seek static profiles x 3 = x 3 (a) that satisfy the equations of 
motion from ( |A.2|) . We rescale the variables 

r = r y, x 3 = — z , (A. 3) 

and we introduce the notational short hands 

H = y 4 - cosh 2 r] , D=y A -l. (A.4) 
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The world-sheet action S takes the form [z! = d y z 



S = K J\y^+Dz», (A.5) 

with K = y/XTT, where T is the extension of the Wilson loop in the t-direction. The Euler 
Lagrange equations of motion imply |^ = q = const., which leads to 

z' 2 = q 2 , / ,/- C ° Sh2V 2 . (A.6) 

We now classify the solutions to these equations of motion which begin from y = A. 

A.l Solutions in the A > \J cosh rj regime 

1. To have a solution stretching between y = A and y — 1, the only allowed value for q is 

q 2 = sinh 2 77 , (A. 7) 

which leads to 



sinh//— -. (A. 8) 



1 

Integration of this equation gives the drag solution 

z(y) = const. — sinh?] [arctan(y) + arccoth(?/)] . (A. 9) 

This is the solution in the rest frame of the quark and its string. In the rest frame of the 
medium, (|A.9p is multiplied by the Lorentz contraction factor yl — v 2 = 1/ coshr/, and agrees 
with the solution of Ref. |23|, 24]. For the solution ( |A.8|) , the action (|A.5|) takes the form 

dy. (A.10) 



S = VXTT J 



Since T is the proper time in the rest frame of the quark, it is related to the "laboratory" 
time in the rest frame of the medium via T = t^/ cosh?]. Hence, ( |A.10| ) is the relativistic 
boost of the action of a static quark. This makes it the natural choice for the subtraction 
term in our analysis of the quark- ant iquark potential in Sections and 5. We saw in Fig. |2] 
that, with this choice of subtraction, the quark- ant iquark potential at small L is independent 
of the velocity of the medium, as is desirable on physical grounds. 

2. Solutions of (|A.6|) , which stop at y\ = a/ cosh 77, exist for values of q satisfying 

< q 2 < sinh 2 77 . (A. 11) 

The actions for such solutions are 



So = K f iy Jtf^. (A . 12) 

J vi V V 1 ~ ? 

These solutions describe "mesons" made from a heavy quark and a light quark with the light 
quark dragging behind the heavy quark. 



- 50 - 



A. 2 Solutions in the ^/cosh rj > A regime 



All solutions in this regime stretch between y = A and the horizon y 
cosh 2 1] < 0, since y < A. Requiring z' 2 > 0, one finds from ( |A.6|) that 



1. In this regime, y A — 



y A - (1 + q 2 ) 

This condition can be realized in two different ways: 

1. Solutions with time-like world sheet. 
Eq. ( |A.13| ) can be satisfied for 



< 



q 2 > 



y 



[l + q 2 ) <0 



The action for these solutions of ( |A.6| ) is time-like 



l + q > A 



S = K j dy 



' cosh r] — y 4 
1 + q 2 — y 4 



(A.13) 



(A. 14) 



(A.15) 



For the value q 2 = sinh r], this action coincides with (|A.10| ) 



2. Solutions with space-like world sheet. 
Eq. ( |A.13| ) is also satisfied for 



q 2 < , y 4 - (1 + q 2 ) > ► 1 + q 2 < A 4 . 

In this case, q = ip is purely imaginary. The equation of motion ( |A.6|) becomes 



p a/ cosh 2 rj — y A 

y 4 -i ^/yT^rrr 



(A.16) 



(A.17) 



which has well-defined solutions for real values p > 0. The action is that of a space-like world 
sheet and is imaginary 



S = iKj\yl™ Sh2 1-»' 



y 4 +p 2 -l 



(A.18) 



In the calculation in Section [3.2| of the expectation value of the light-like Wilson loop that defines 
the jet quenching parameter, we used ( [A. 18| ) with p = as the L-independent subtraction term 
because it satisfies ( |3.43|) . 
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B. Time-like Wilson loop with dipole parallel to the wind 

In Section |], we calculated the static gg-potential for (almost) all dipole orientations with respect 



to the wind. While the parameterization (5T) of the two-dimensional world sheet used there 



is applicable for arbitrarily small angles 9, it is not applicable for 9 = 0, the case where the 
dipole is parallel to the wind. In this appendix, we repeat the calculation of Section |5] with a new 
parametrization which works for 9 = and in fact for < 9 < n/2 but which does not work for 
9 = tt/2. 

We start again from the boosted metric ( p. 13 ), but in contrast to (|5.1| ), we parametrize the 
world sheet by 

r = t, a = x 3 , x 2 = const, Xi=X\{a), r = r(a) . (B.l) 

The role of Xi and x% are interchanged in this parametrization relative to that in Section |^. 
We define dimensionless coordinates 

* = -, * = - = ^ * = (B.2) 

and drop the tilde. The boundary conditions on y(cr) and w(a) then become 

y (± l -cos6 J = A, w ( ±-cos9 J = ±-sin0. (B.3) 



2 J V 2 / 2 

-1/2 



The Nambu-Goto action takes the form S(C) = y/XTT fJ daC, with the Lagrangian now given by 



C=^J(y*-l) + (y* - cosh 2 ri)w>* + V - V y* , (B.4) 

where y' and w' denote derivatives with respect to a. The constants of the motion are 

H = £-y'— -w'— - ^ - 1 = q 
dy' dw' C 

dC y 4 - cosh 2 r/ , 

d^ = c z = p - (R5) 

Note that the constants of the motion p and q here are not the same as in Section ^. The equations 
of motion can be written in the form 

q 2 y' 2 = ( 4 V " " \ X [(y* - cosh 2 r,)(y 4 - 1 - g 2 ) - p\y* - 1)] , (B.6) 
\?/ 4 — cosh 7] J 

q V= P 2 { /~] 2 (B.7) 
— cosh rj J 
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Since y' becomes singular at y 2 = cosh 77, the turning point y c which satisfies 



(yt - cosh 2 V )(y 4 - I - q 2 ) - V \v\ ~ 1) = (B.8) 
must also fulfill the condition 

\J cosh r] < y c < A. (B.9) 
The constants q and p are related to the values of I and 6 via 

1 ■ ft [ Adw j 
- sin 6 = / -j—dy 



2 L dy 



oo 



/'/ ^ =, (B.10) 

(y 4 - cosh 2 77) (y 4 - 1 - g 2 ) - p 2 (t/ 4 - 1) 



II r 



1 a f°° y ~ cosh ^ d V 
cos6 = q I - ; (B.ll; 

(y 4 - cosh 2 rj) (t/ 4 - 1 - q 2 ) - p 2 {y 4 - 1) 



2 V* l/ 4 -l 

We now see that the expressions ( |B.8|) , ( |B.10| ) and (|B.12|) differ from their analogues in Section 5 
simply by exchanging q and p. 

Results with < 6 < tt/2 can be obtained with either the parametrization of the string world 
sheet in this Appendix or that in Section 5. Let us now specialize to 9 = 0, the case that cannot 
be handled with the parametrization of Section 5. We see from ( [B.10| ) that 6 = corresponds to 
p = 0, which means that the turning point which solves (p.8|) is given simply by 



y 4 c = l + q 2 . (B.12) 

Now, the condition ( |B.9| ) becomes a restriction on the allowed values of q: 

sinh 2 r] < q 2 < A 4 - 1 . (B.13) 

In contrast to the situation for any nonzero value of 9, when 6 = the constant q cannot be taken 
to zero: if we were to choose q < sinh 77, then y 4, would hit cosh 2 ?], at which point y' — > 00 and 
below which y' is imaginary. Instead, when we choose q > sinh?] we find a solution in which y 
reaches a turning point at y c and safely begins to ascend, never reaching these pathologies. 

Notice that there is "almost" another choice of y c : y c = a/ cosh 77 does satisfy ( p.8| ) with p = 0, 
but it just barely fails to satisfy ( |B.9| ). For arbitrarily small but nonzero values of p, however, if 



q < sinh?? there is a legitimate turning point at a y c just above ycosh??, and a solution can be 
found. So, it is only for 6 = and hence p = that there is an inaccessible range of small values 
of q. If we sit at a q which is less than sinhr? and take p — > 0, what happens to the solution is 
that the string world sheet develops a cusp at its turning point. For p small but nonzero, the shape 
of the function y(<r) near its minimum looks like a very slightly rounded "V", with the amount of 
rounding controlled by p. So, for p = there is no solution in this regime of small q. We shall see 
momentarily that this regime of q corresponds to a part of the unstable higher energy branch of 
solutions. 
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Let us return to the case with 9 = p = keeping q > smhr/. The solution evidently has 
w' = throughout. Furthermore, although there is no reason of symmetry to expect it, given the 
wind blowing in the £3 direction, we can see that the solution y(cr) will be a — > —o symmetric. 
This follows from the fact that y' 2 depends only on y, not explicitly on a, and from the symmetric 
boundary condition (|B.3|) . This means that the descending half of the y(o) curve and the ascending 
half must have the same shape, implying that the turning point at which y = y c must be at a = 0, 



half way between the boundaries at which the boundary condition (|B.3|) fixes y. q can be determined 
in terms of I from ( B.ll ), which with p = is simply the equation | = J 2 da and becomes 



i = 2,/"V V/ ^ 13 ™* 5 



with y c given by ( B.12| ). The action can be written as 

\/ y 4 — cosh 2 rj 



[B.U) 



3(1) = VXTT ! dy 



[B.15) 
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from which the quark- ant iquark potential E(L) can be obtained as in (|3.27|) using the same sub- 
traction So given in ( |3.28| ). In plotting Fig. we have used ( |B.14j ) and ( |B.15| ) to evaluate I and 
E for q > sinh?]. We see that the stable, lower energy branch of solutions is similar to those we 
have obtained previously. For this branch of solutions, the 9 — > limit is smooth. The inaccessible 
range of q, namely q < sinh rj, where as described above the string world sheet develops a cusp in 
the 9 — > limit, corresponds to the "missing parts" of the unstable high energy branch of solutions 
in Fig. R. 
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